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Non- homogenous superconductors and non- homogenous superfluids appear in a variety 
of contexts which include quark matter at extreme densities, fermionie systems of cold 
atoms, type-II cuprates and organic superconductors. In the present review we shall 
focus on the properties of quark matter at high baryonic density which can exist in the 
interior of compact stars. The conditions that are realized in this stellar objects tend 
to disfavor standard symmetric BCS pairing and may be in favor of a non-homogenous 
color superconducting phase. We discuss in details the properties of non-homogenous 
color superconductors and in particular of crystalline color superconductors. We 
also review the possible astrophysical signatures associated with the presence of 
non-homogenous color superconducting phases within the core of compact stars. 
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I. INTRODUCTION 



Ideas aboii t color supercondiicting (CSC) matter d ate back to more than 30 years ago (jBailin and Lovd . 11984 ; 



Barroisl 1977: Collins and Perrv, 1975; 


Frautschi, 1978). but only recentb 


of consideration (for recent reviews see 


Alford (2001):lAlford et all ( 


2008^ 


RaiaeoDal and Wilczek ( 


20001): Rischkd (20041): andlSchafeil (2003b 


)). c 



analogous of the standard electromagnetic superconductivity and is believed to be the ground state of hadronic matter 
at sufficiently large baryonic densities. At very high density the naive expectation, due to asymptotic freedom, is that 



quarks f orm a Fermi sp here of almost free fermions. However, Bardeen, Cooper and Schrieffer (BCS) (jBardeen et al. 
1957313; [Cooce3,[ 19561 ) have shown that the Fermi surfaces of free fermions are unstable in presence of an attractive. 



arbitrary small, interaction between fermions. I n Quantum-Chromodynamic s (QCD) the attractive interaction be- 
tween quarks can be due to instanton exchange (jSchafer and Shurvakl . 119981 ). at intermediate densities, or to gluon 



3 



exchange in the 3 color channel, at higher densities. Therefore, one expects that at high densities quarks form a 
coherent state of Cooper pairs. 

It shou l d be n oticed that the mentioned old papers ( Bailin and Love . 1984 : Barroij 119771: ICollins and PerrvL [lOTSi: 
Frautschil Il978l ) were based on the existence of the attractive 3 color channel and on analogies with ordinary super- 
conductors. The main result of these analyses w as th at quarks form Coo per pairs with a gap of order a few MeV. In 
more recent times two papers ( Rapp et al . 19981 ) and ( Alford et al . 1998[) . have brought this result to question. These 
authors considered diquark condensation arising from instanton-mediated interactions and although their approxima- 
tions are not under r igorous qua ntitative control, the result was that gaps can be as large as 100 MeV. However, as 
point ed out by Hsu Jhs^ . l200n[ ): " this is an unfair characterization of Bailin and Love 's results ^Bailin and Lov^ . 
198 A ). A v a lue of the strong coupling large enough to justify the instanton liquid picture of lAlford et al. 1(1994 ) and 
Raw et al.l H994 } also yields a large gap when substituted in Bailin and Love 's results. After all, instantons are sup- 
pressed by an exponential factor exp(— 27r/as). Bailin and Love merely suffered from the good taste not to extrapolate 
their results to large values of Ug!^^ 

Color superconductivity offers a clue to the behavior of strong interactions at very high baryonic densities, an 
issue of paramount relevance for the understanding of the physics of compact stars and of heavy ion collisions. In 
the asymptotic regime it is possible to understand the structure of the quark condensate from basic considerations. 
Consider the matrix element 



(oiCV-fjo), 



(1) 



where tpf^, V'jt represent the quark fields, and a, /? = 1, 2, 3, s, t = 1, 2, i, = 1, • • • , Nf are color, spin and flavor indices, 
respectively. For sufficiently high quark chemical potential, /i, color, spin and flavor structure can be completely fixed 
by the following arguments: 

• Antisymmetry in color indices (a, /?) in order to have attraction. 

• Antisymmetry in spin indices (s,t) in order to have a spin zero condensate. The isotropic structure of the spin 
zero condensate is favored with respect to the spin 1 (or higher spin) condensate since a larger portion of the 
phase space around the Fermi surface is available for pairing. 

• Given the structure in color and spin, Fermi statistics requires antisymmetry in flavor indices. 

Since the quark spin and momenta in the pair are opposite, it follows that the left(right)-handcd quarks can pair 
only with left(right)-handed quarks. Considering three-flavor quark matter at large baryonic density, the so-called 
color-flavor locked (CFL) phase turns out to be thermodynamically favored, with condensate 



|0) 



-(0|^r«<i^|0)«AcFLVe"^%/ 



1=1 



(2) 



where Acfl is the pairing gap, we have suppressed spinorial indices and neglected pairing in the color sextet channel. 
Pairing in the color sextet channel is automatically induced by the quark color structure, but the condensate in 
this channel is much smaller than in th e color anti-t r iplet c hannel and in most cases can be neglected. The CFL 
condensate in Eq. ([2|) was introduced in lAlford et al\ (Il999cl ) and the reason for its name is that only simultaneous 
transformations in color and in flavor spaces leave the condensate invariant. The symmetry breaking pattern turns 
out to be 



S'C/(3)c ® 5'C/(3)l ® SU{3)b^ U{1) 



SU{3) 



c+L- 



-R ® Z2 ; 



(3) 



where SU{3)c+l+r is the diagonal global subgroup of the three SU{3) groups and the Z2 group means that the 
quark fields can still be multiplied by -1. According with the symmetry breaking pattern, the 17 generators of chiral 
symmetry, color symmetry and U{1)b symmetry are spontaneously broken. The 8 broken generators of the color 
gauge group correspond to the 8 longitudinal degrees of freedom of the gluons and according with the Higgs- Anderson 
mechanism these gauge bosons acquire a Meissner mass. The diquark condensation induces a Major ana-like mass 
term in the fermionic sector which is not diagonal in color and flavor indices. Thus, the fermionic excitations consist 
of gapped modes with mass proportional to AcflIII The low-energy spectrum consists of 9 Nambu-Goldstonc bosons 



^ This is a feature of all homogenous superconducting phases: the fermionic excitations which are charged with respect to the condensate 
acquire a Majorana-like mass term proportional to the pairing gap. 
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(NGB) organized in an octet, associated with the breaking of the flavor group, and in a singlet, associated with the 
breaking of the baryonic number. For non- vanishing quark masses the octet of NGBs becomes massive, but the singlet 
NGB is protected by the U{1)b symmetry; it remains massless and determines the sup erfluid properties of the CFL 
phase. The effective theory des c ribing the NGBs for the CFL phase has been studied in (jCasalbuoni and Gattd . ll999l: 



Son . 2002 : Son and Stephanovl 2000a). The CFL condensate also breaks the axial U{1)a symmetry; given that at 



very high densities the explicit axial symmetry breaking is weak, one has to include the corresponding pseudo-NGB 
in the low-energy spectrum. 

After the first attempts with instanton-induccd interaction many authors tried various approaches in order to 
calculate the magnitude of gap param eters in the CSC phases (for references see the review by Rajagopal and 
Wilczek (jRaiag opal and Wilczekl . l2000f ) ) . Dealing with QCD the ideal situation would be if these kind of calculations 
could fall within the scope of lattice gauge theories. Unfortunately, lattice methods rely on Monte Carlo sampling 
techniques that are unfeasible at finite density because the fermion determinant becomes complex. Alt hough various 



appro ximation schemes have b een developed, fo r insta nce, Taylor expansion in the chemical potential (|Allton et al. 



2003 ). reweighting technique s ( Fodor and Katz . 200^ . analytical continua tion of calculation employ ing imaginary 



baryo nic chemical potential ( Alford et all 1999b : Roberge and Weissl . Il986i) or heavy Wilson quarks ( Fromm et al 



20121 ). no definite results have been obtained so far for large values of the baryonic chemical potential and physical 



quark masses. 

In the absence of suitable lattice methods, quantitative analyses of color superconductivity have followed two distinct 
paths. The first path is semi-phenomenological, and based on simplified models. The main feature of these models is 
that they should incorporate the most important physical effects but being at the same time tractable within present 
mathematical techniques. All these models have free parameters that are adjusted in such a way to give rise to a 
reasonable vacuum physics. 

Examples of these kind of techniques include Nambu-Jona Lasinio (NJL) models in which the interaction be- 



tween quarks is replaced by a four-fermion in teraction originating from instanton exchange (jAlford et ali . 11998 : 



Berges and Raiagopal . 1999 : Raw et aZ.I. 19981) or where the four-fermi interaction is mode led by that induced by 



single- gluon exchange (jAlford et all Il999al 119981 ) . Random m atrix models have been studie d inlVanderheyden and Jackson 
(|2000[ ) and instanton liquid models have been investigate d in ICarter and Diakonovl (fl999[) : iRapp et al\ (l200dl200lh . 
while renormalization group methods have been used in lEvans et all (jl999 ) and Schafer and Wilczek (1999b). Al- 
though none of these methods has a firm theoretical basis, all of them yield results all in fairly qualitative agreement. 
This is probably due to the fact that what really matters is the existence of an attractive interaction between quarks 
and that the parameters of the various models arc chosen in such a way to reproduce the chirally broken ground 
state. The gap parameter evaluated within these models varies between tens of MeV up to 100 MeV. The critical 
temperature is typically the same found in normal superconductivity, that is about one half of the gap. 

The second path starts from first principles and relies on the property of asymptotic freedom of QCD. Various 
results have been obtained star ting from the QCD action, emp l oying renormalization group techniques and through 
the Schwinger-Dyson equation ([Brown et all \200d lEvans et al\.\200d iHong et all l2000l : iPisarski and Rischkd . l2000at 
Schafer and Wilczekl . Il999cl: ISonl . ll999t ). In particular. Son (jSonl Il999l ) , using the renorr nalization group near the Fermi 



surface has obtained the asymptotic form of the gap. However, it has been argued in Raiagopal and Shuster ( 2000l ) 
that the weak-coupling calculations employed in this case are valid only for extremely large densities corresponding 
to chemical potentials larger than 10^ MeV. Moreover, also in this second path, it is not possible to control the 
approximation, because we do not know how to evaluate higher order corrections. Interestingly enough, extrapolating 
the results from asymptotic densities to values of quark chemical potential of order (400 — 500) MeV, one obtains a 
magnitude of the gap in agreement with the results of more phenomenological approaches. 

The result of the analyses of the above-mentioned methods is that the CFL phase is the thermodynamically favored 
state of matter at asymptotic densities. Qualitatively one can understand this result considering that in the CFL phase 
quarks of all thrce-fiavors participate coherently in pairing. Since superconductivity is a cooperative phenomenon, the 
larger the number of fcrmions that participate in pairing, more energetically favored is the superconducting phase. 

In the description of color superconductivity one has to deal with various scales, the chemical potential the gap 
parameter, which wc shall gcnerically indicate with A, the strange quark mass Ms and the screening/damping scale 
g^, where g is the QCD coupling. One typically has that fi ^ gfi ^ whereas the strange quark mass can be 
considered as a free parameter, although in some models it can be computed self-consistently. 

QCD at high density is conveniently studied through a hierarchy of effective field theories, schematically depicted 
in Fig. [TJ The starting point is the fundamental QCD Lagrangian, then one can obtain the low energy effective 
Lagrangian thro ugh different methods. One way is to integrate out high-energy degrees of freedom as shown in 



Polchinskil (jl992t ). The physics is particularly simple for energies close to the Fermi energy where all the interactions 
are irrelevant except for a four-fermi interaction coupling pair of fermions with opposite momenta. This is nothing 
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P+5 



P+A 



-'QCD 



■'NGB . 



FIG. 1 Schematic representation of the hierarchy of effective Lagrangians characteristic of high density QCD. 



but the i nteraction Rivhig; rise to BCS condensation, which can be d e scribed using the HiRh Density Effective Theory 
(HDET) (jBeane et aLl . l2000HCasalbuoni et ad [200 lbl:lHond . l2000allbl : lNa"rdulli 120021: ISchafeill2003al) . see Sec. Illl.B.ll 
for a brief description of the method. The HDET is based on the fact that at vanishing temperature and large 
chemical potentials antiparticlc fields decouple and the only relevant fermionic degrees of freedom are quasi-particles 
and quasi-holes close to the Fermi surface. In the HDET Lagrangian the great advantage is that the effective fermionic 
fields have no spin structure and therefore the theory is particularly simple to handle. 

This description is supposed to hold up to a cutoff + J, with 6 smaller than the Fermi momentum, ^ but 
bigger than the gap parameter, i.e. A ^ (5 ^ P^ . Considering momenta much smaller than A all the gapped 
particles decouple and one is left with the low energy modes as NGBs, ungapped fermions and holes and massless 
gauged fields according with the s ymmetry breaking schern e . In the case of CFL arid othe r CSC phases, such e f Tectiv e 
Lagrangians have been derived in ICasalbuoni et al\ (|2000t) : ICasalbuoni and Gattol (|1999[) : and iRischke et al\ (|200l[ ). 
The parameters of the effective Lagrangian can be evaluated at each step of the hierarchy by matching the Green's 
functions with the ones evaluated at the upper level. For the CFL phase the effective Lagrangian of th e superflui d 
mode associated with the breaking of C/(1)b may also be determined by symmetry arguments alone as in lSonI (20021). 

In the high density limit one neglects the quark masses, but in realistic situations the quark chemical potential 
may be of the same order of magnitude of the strange quark mass. The typical effect of quark masses is to produce a 
mismatch between Fermi surfaces. Neglecting light quark masses and assuming (for simplicity) that quarks have all 
the same chemical potentials, the Fermi spheres have now different radii 



11^ 



AP P^ 



(4) 



where Mg is the constituent strange quark mass. Thus, increasing for a fixed value of fi, increases the mismatch 
between the Fermi surface of strange quarks and the Fermi surfaces of up and down quarks (which in this simplistic 
case are equal). 

The standard BCS mechanism assumes that the Fermi momenta of the fermionic species that form Cooper pairs 
are equal. When there is a mismatch it is not guaranteed that BCS pairing takes place, because the condensation 
of fermions with different Fermi mome nta has a free-energy cost. For small mismatches there is still condensation 
(jChandrasekharl . Il962l : IClogstonl Il962l ). and in the case at hand it means that the CFL phase is favored. However, 
for large values of the strange quark mass the assumptions leading to prove that the favored phase is CFL should be 
reconsidered. According with Eq. (|4]) if the strange quark mass is about the quark chemical potential, then strange 
quarks decouple, and the corresponding favored condensate should consist of only up and down quarks. With only 
two flavors of quarks, and due to the antisymmetry in color, the condensate must necessarily choose a direction in 
color space and one possible pairing pattern is 



(OlCL^filO) (X A2sce"^%-3, a,peSUc{2) i,jeSU{2)i 



(5) 



This phase of matter is known as two-flavor color superconductor (2SC) and A2SC is the corresponding gap param- 
eter. This phase is characterized by the presence of 2 ungapped quarks, qub,Qdb and 4 gapped quasi-particles given 



by the combinations 



que 



and 



qdg of the quark fields, where the color indices of the fundamental repre- 



sentation 1,2,3 have been identified with r,g,b (red, green and blue). In case massive strange quarks are present 
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the corresponding phase is na med 2SC+S and eventually strange quarks may by themselves form a spin-1 condensate 
(IPisarski and Rischkd . boOObI ). 

In the 2SC phase the symmetry breaking pattern is completely different from the three-flavor case and it turns out 
to be 



SU{3)c (E) SU{2)l'E) SU{2)r ® U{1)b ^ S'f/(2), ® SU{2)l ® SU{2)r ® ?7(1)b ® ^2 



(6) 



The chiral group remains unbroken, meaning that there are no NGBs. The original color symmetry group is broken to 
SU{2)c and since three color generators are unbroken, only five gluons acquire a Meissner mass. Even though U{\)b 
is spontaneously broken there is an unbroken U{l)g global symmetry, where B is given by a combination of B and 
of the eighth color generator, Ts, playing the same role of the original baryonic number symmetry. In particular, this 
means that unlike CFL matter, 2SC matter is not superfiuid. One can construct an effective theory to describe the 
emergence of the unbroken subgroup SU(2)c and the low energy excitations, much in the same wa y as one builds up 
chiral effec tive Lagrangian with effective fields at zero density. This development can be found in ICasalbuoni et al 
(|2000l ) and lRischke et all (|200ll ). 



The problem of condensation for imbalanced Fermi momenta is quite general and can be des cribed assuming that 



different populations have different chemical potentials, see e.g. lAlford and Raiagopall (|2002h . We shall present a 
discussion of a simple two-level system in Sec. III.Al which allows to clarify the main effects of a chemical potential 
difference, Sji. 

Regarding quark matter a similar behavior can be found when one considers realistic condition, i.e. conditions 
that can be realized in a compact stellar object (CSO). This is a real possibility since the central densities for these 
stars could reach 10^^ g/cm^, whereas the temperature is of the order of tens of keV, much less than the critical 
temperature for color superconductivity. Matter inside a CSO should be electrically neutral and in a color singlet 
state. Also conditions for /3-equilibrium should be fulfilled. As far as c olor is concerned, it is possible to impose a 
simpler condition, that is color neutrality, since in lAmore et ali (|2002l ) it has been shown (in the two-flavor case) 



that there is a small free-energy cost in projecting color singlet states out of color neutral ones. If electrons are 
present (as generally required by electrical neutrality) the /^-equilibrium condition forces the chemical potentials of 
quar ks with different electric char ges to be different, thus Eqs. (j4]) are no more valid and more complicated relations 
hold (lAlford and Raiagopall [200^ . The effect of the mass of the strange quark, /3-equilibrium and color and electric 
neutrality, is to pull apart the Fermi spheres of up, down and strange quarks and many different phases may be 
realized, depending on the parameters of the system. Besides the above-mentioned standard 2SC and 2SC-I-S phases, 
the two-fl a vor su perco nducting phas e 2SCus , with pairing between up and strange quarks can be favored; see e.g. 
lida et ali (|200-i ) and iRuester et ali (|2006al ) for different pairing patterns. For very lar ge mismatche s among the 



lida et ali l|zUU4l l and IKuester et ali l|zUUbal i tor ditterent pairmg patterns, l^or very lar ge mismatcne s among the 
three flavors of quark s only the inte r -spec i es single-flavor spin-1 pairing may take place (lAlford et al. , 20031 1998t 



i P ^ 

BaiUn and Lwi Il979l: iBuballa et d\ . l2003l : ISchafeil . I2OOOI: ISchmittl . 
(|2008l ) for an extended discussion on these topics 



2005HSchmitt et ad . 12003) . see e.g. lAlford et al. 



Sec 



In the 2SC phase, the above conditions tend to separate the Fermi spheres of up and down qua rks and, as discussed in 
I.B[ for |i5/i| = A gapless mo des appear. The corresponding phase has been named g2SC (|lluang and Shovkovv . 
Shovkovv and Huangl . l2003l ) . with "g" standing for gapless. The g2SC phase with pairing between up and down 



2003 



quarks has the same condensate of the 2SC phase reported in Eq. (O, and therefore the ground states of the 2SC and 
of the g2SC phases share the same symmetry. However, these two phases have a different low energy spectrum, due 
to the fact that in the g2SC phase only two fermionic modes are gapped. The g2SC phase is energetically favored 
with respect to the 2SC phase and unpaired quark matter in a certain ra nge of values of the four-fer mi interaction 
strength when one considers /3-equilibrium, color and electrical neutrality ([Shovkovv and Huangl . l2003i ). 

Pinning down the correct ground state of neutral quark matter in /3-equilibrium is not simple because another 
difficulty emerges. This pro blem, which is already p resent in the simple two-level system discussed in Sec. III.Al has 
a rather general character ( Afford and Wan^ . 20051) . and is due to an instability connected to the Meissner mass. 
In particular, when \Sfj,\ = A the system becomes magnetically unstable, meaning that the Meissner mass becomes 
imaginary. In the 2SC phase the color group is broken to SU{2)c and 5 out of 8 gluons acquire a mass. Four of 
these masses turn out to be imaginary in the 2SC phase for A/^/2 < Sfi < A, thus in this range of the 2SC 
phase is chromo-magnetically unstable ( Huang and Shovkovvl [2004al lb[). Increasing the chemical potential difference 
the instability gets worse, because at the phase transition from the 2SC phase to the g2SC phase all the five gluon 
masses become pure imaginary. 



An analogous phenomenon arises in three-flavor quark matter because the gapless CFL fg CFL) phase (lAlforc 



Fukushimal . 
between up 



20051 ) turns out to be chromo-magnetically unstable (|Casalbuoni et al 



et al. 



2005bl : 



2005bL 12004 12005d iFukushima et al 

20051) . The gCFL phase has been proposed as the favored ground state for sufficiently large mismatch 



down and strange quarks and occurs in color and electrically neutral quark matter in /3-cquilibrium 
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for > A. Some properties of the gCFL phase and a brief discussion of the corresponding instabiUty will be 

presented in Sec. IIII.AI 

There is a variety of solutions that have been proposed for the chromo-magnetic instability and that can be realized 
depending on the particular conditions considered. As we have already discussed, the chromo-magnetic instability is 
a serious problem no t only for the gapless phases (g2SC and gCFL) but also for the 2SC phase. In the latter case, 
it has been shown in lGorbar et al\ (|2006al ) that vector condensates of gluons with a value of about 10 MeV can cure 
the instability. The corresponding phase has been named gluonic phase and is characterized by the non-vanishing 
values of the chromo- electric condensates (A"^) and (A^) which spontaneously break the S0(3) rotational symmetry. 
It is not clear whether the same method can be extended to the gapless phases. The chromo-magnetic instability 
of the gapped 2SC phase can also be removed by the formation of an inhomogeneous condensate of charged gluons 
(jFerrer and de la Inceral . 120071 ) . 

For the cases in which the chromo-magnetic instability is related to the presence of gapless modes, in iHond (l2005l) 
the possibility that a secondary gap opens at the Fermi surface is studied. The solution of the instability is due to a 
mechanism that stabilizes the system preventing the appearance of gapless modes. However, the secondary gap turns 
out to be extremely sma ll and at temperatures typical of CSOs it is not able to fix the chromo-magnetic instability 
(lAlford and Wang| .' l2006l) . 

The imaginary value of the Meissner mass can be understood as a tende ncy of the system toward a non- homogeneous 
phase (jGubankova et all l2010l : iHond . 120051: llida and Fukushimal . l2006l ). This can be easily seen in the toy-model 
system discussed in Sec. III. Al for the case of a U{1) symmetry, where one can show that the coefficient of the gradient 
term of the lo w energy fluctuations ar ound the ground state of the effective action is proportional to the Meissner 



mass squared ( Gubankova et al . 2010[ ) 



For three-flavor quark matter two n on-homogenous super c onducting phases have been pro posed. If kaon conden- 
sation takes place in the CFL phase (jBedaaue and Schaferl . l2002t iKaplan and Reddvl . 120021) , the chromo-magnetic 
instability might drive the system toward a non-homogeneous state where a kaon condensate current is generated, 
balanced by a counter-propagating current in the opposi te direction carr ied by gapless quark quasi-particles. This 
phase of matter, named curCFL-_ftr*', has been studied in iKrvievskil (2008|) and turns out to be chromo-magnetically 
stable. 



Th e second possibility is the crystalline color superconducting (CCSC) phas e ( Alford et clI\. 2001 Bowers et al 



2001 



Bowers and Rajagopaj , 



2002! iLeibovich al\ . (200 1[ 



2002tlCasalbuoni et adf200ll2002b[|2005ail2001all2002cll2003HKundu and Raiagopal . 



Mannarelli et all l2006a l , which is t he QCD analogue of a form of non-BCS pairing 



first proposed by Larkin, Ovchinnikov, Fulde and Ferrell (LOFF) ( Fulde and Ferrell . 1964t Larkin and Ovchinnikov 
19641 ). This phase is chromo-magnetically stable as we shall discuss in Sec. III.Fl for two-fiavor quark matter, and in 
Sec. nil. B. 41 for three-fiavor quark matter. The condensate characteristic of this phase is given by 



1=1 



(7) 



which is similar to the condensate reported in Eq. ^ but now there are three gap parameters, each having a 
periodic modulation in space. The modulation of the /'th condensate is defined by the vectors q™, where m is the 
index which identifies the elements of the set {qj}. In position space, this corresponds to condensates that vary 
like X)in *^^P(2iQ'" ' ''')■, meaning that the <7""s are the reciprocal vectors which define the crystal structure of the 
condensate. 

In Sec, nil we sh all discuss various properties of the the two- flavor CCSC phase. This phase, first proposed in 
Alford et al\ (|200lh . corresponds to the case where only one gap parameter is non- vanishing. The chromo-magnetic 
stability of a simple two-flavor periodic structure with a gap par ameter modulated by a single plane wave (here- 
after we shall refer to t his phase as Fulde-Ferrell (FF) structure I Fulde and FerreUl . 19641 )) has been considered in 
Giannakis et al. ( 2005h : Giannakis and RenI (j2005al lbl) where it has been shown that 



• The presence of the chr omo-magnetic in stability in g2SC is exactly what one needs in order that the FF phase 
is energetically favored (jGiannakis and Re n. 2005al) . 



The FF phase in the two-fl avor case has no chromo-magnetic instability (thou gh it has gapless modes) at least 
in the weak coupling limit (jGiannakis et al\ . l2005l : I Giannakis and "RCTl . l2005bl) . 



The stability of the FF phase in the strong coupling case has been studied in iGorbar et aO (|2006bh . where it is 
shown that for large valu es of the gap parameter the FF phase cannot cure the chromo-magnetic instability. In 
Nickel and Buballa (I2OO9I) it has been questioned whether among the possible one-dimensional modulations, the 
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periodic LOFF solution is the favored one. As we shall discuss in Sec. III.Gi it is found that for two-flavor quark 
matter, a solitonic-like ground state is favored with respect to FF in the range of values 0.7A ^ Sfi < 0.78 A. 
However, at least in weak coupling, the FF phase is not the crystalline structure one should compare to. The FF 
phase is slightly energetically favored with respect to unpaired quark matter and 2SC quark matter for < Sfj, < 

0.7 54 A, but more complicated c rystalline structures have larger condensation energies in a larger range of values of 
5fi ( Bowers and Raiagopall [2002h . 



The stability analysis of the three-flavor CCSC phase is discussed in Sec. IIII.B.41 where we repo rt on the results 



obtai ned for a simple structure made of two plane waves by a Ginzburg-Landau (GL) expansion (jCiminale et at. 



200^ . This particular three-flavor CCSC phase turns out to be chromo-magnetically stable, but the stability of more 



complicated crystalline structures has not been studied, although by general arguments they are expected to be stable, 
at least in the weak coupling limit. 

Whether or not the crystalline color superconducting phase is the correct ground state for quark systems with 
mismatched Fermi surfaces has not yet been proven. In any case it represents an appealing candidate because in 
this phase quark pairing has no energy cost proportional to i5/i. The reason is that pairing occurs between quarks 
living on their own Fermi surfaces. However, this kind of pairing can take place only if Cooper pairs have nonzero 
total momentum 2g and therefore it has an energy cost corresponding to the kinetic energy needed for the creation of 
quark currents. Moreover, pairing can take place only in restricted phase space regions, meaning that the condensation 
energy is smaller than in the homogeneous phase. The vector q has a magnitude proportional to the chemical potential 
splitting between Fermi surfaces, whereas its direction is spontaneously chosen by the system. In case one considers 
structures c omposed by a set of v e ctors |qr|, one has to find the arr angement that minimizes the free-energy of 



the system ( Bowers and Raiagopal . 2002t Raiagopal and Sharmal . [2006h . This is a rather complicated task which is 



achieved by analyzing some ansatz structures and comparing the corresponding free-energies. We shall report on this 
subject in Sec. III. CI and III. El for the two-flavor case and in Sec. Illl.Cl for the three-flavor case. 

Summarizing, we can say that the state of matter at asymptotic densities is well defined and should correspond 
to the CFL condensate. At intermediate and more realistic densities it is not clear which is the ground state of 
matter. Our knowledge of the phases of matter can be represented in the so-called QCD phase diagram, which is 
schematically depicted in Fig. [51 At low-density and low temperature quarks are confined in hadrons but increasing 
the energy scale quarks and gluons degrees of freedom are liberated. At high temperature this leads to the formation 
of a plasma of quarks and gluons (QGP), while at large densities matter should be in a color superconducting phase. 
Apart from the phases we have discussed other possibilities may be realized in the density regime relevant for CSO. 
Here we only mention that rec ently it has been proposed one more candidate phase, the so-called quarkyonic phase 
( McLerran and Pisarski . 2007t ). which is characterized by a non- vanishing baryon number density and found to be a 



candidate phase at least for a large number of colors. Another possibility is that the constituent value of the strange 
quark mass is so small that the CFL phase is the dominant one down to the phase transition to the hadronic phase. In 
this case, a rather interesting possibility is that there is no phase transition betw een the CFL phase and the h adronic 



phase (hyper- nuclear matter) . in the so-called quark-hadron continuity scenario (jSchafer and Wilczekl . Il999a^ . 

In Section IV we shall discuss whether the presence of a non-homogenous color superconducting phase within the 
core of a compact star may lead to observable effects. The different possible signatures are associated with 

1. Gravitational wave emission. 

2. Anomalies in the rotation frequency (known as glitches). 

3. Cooling processes. 

4. Mass-radius relation. 

Point[T]is discussed in Sec. IIV.AI and relies on the observation that pulsars can be continuous sources of gravitational 
waves if their mass distribution is not axis-symmetric. The large shear modulus characteristic of the CCSC phase 
allows the presence of big deformations of the star, usually called "mountains" . 



Regarding point [H the three-flavor CCSC phase is characterized by an extremely large rigidity (jMannarelli et al 



20071 ). with a shear modulus which is larger than that of a conventional neutron star crust by a factor of 20 to 1000. 
This fact makes the crystalline phases of quark matter unique among all forms of matter proposed as candidates for 
explaining stellar glitches. This topic is detailed in Sec. IIV.BI 

Point [3] is discussed in Sec. IIV.CJ we report the results of a first step toward calculating the cooling rate for neutron 
stars with a CCSC core taken in [Anglani et al. (Hooi), where a simple two plane waves structure was considered. 



Because the crystalline phases leave some quarks at their respective Fermi surfaces unpaired, their neutrino emissivity 
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FIG. 2 (Color online). Schematic phase diagram of strongly interacting matter as a function of the baryonic chemical potential 
and temperature. At low temperatures and low densities matter consists of confined hadrons. At high temperatures quark and 
gluons degrees of freedom are liberated forming the quark-gluon plasma. At low temperatures and very high densities the CFL 
phase is favored. At densities and temperatures relevant for compact stellar objects, the CFL phase may be superseded by 
some different color superconducting phase or by some other phase of matter. The thick blue segment represents the possible 
range of baryonic chemical potential reachable in compact stars. 



and heat capacity are only quantitatively smaller than those of unpaired quark matter (llwamotd . fl980lll98lL not 
parametrically suppressed. This suggests that neutron stars with crystalline quark matter cor es will cool down by the 



direct Urea reactions, i.e. more rapidly than in standard cooling scenarios (jPage et all . l2004f ) 



Point 2] is discussed in Sec. IIV.DI where it is shown that recent observations of very massive compact stars do not 
exclude the possibility that CSO have a CCSC core. 

Finally, studying the damping mechanisms of the radial and non-radial star oscillations (star-seismology) might be 
useful to infer the properties of the CCSC core, and to have information about the width of the various internal layers 
of the star. We are not aware of any paper discussing this topic, which however might be relevant for restricting the 
parameter space of the CCSC phase. 



II. THE TWO-FLAVOR NON-HOMOGENOUS PHASES 



The non-homogeneous two-flavor crystalline color supercon ducting phase is an ext ension to QCD of the phase 
propos ed in condensed matter s ysterr is by Fulde and Ferrell ( Fulde and Ferrelll . 1964 ) and by Larkin and Ovchin- 



nikov (ILarkin and Ovchinnikov 
review 



1964 ) (LOFF). Some aspects of this state have been previously discussed in the 
Casalbuoni and Nardullil (|2004l) . Therefore, we will focus here on recent results, and in particular we discuss 



one of the main properties of these phase, namely its chromomagnetic stability. This important property is not shared 
with other homogeneous gapless color superconducting phases, and therefore it strongly motivates its study. 



A. Mismatched Fermi spheres 

Before discussing the case of two-flavor quark matter, in order to show how gapless superconductivity may arise, it 
is useful to consider the simpler case of a non-relativistic fermionic gas avoiding the formal complications due to flavor 
and color degrees of freedom. We consider a system consisting of two unbalanced populations of different species ipi 
and V'2, with opposite spin, at vanishing temperature, having the hamiltonian density 

H = ^ V^. (-^ - - GM2i^2i^l , (8) 

s=l,2 V / 

where G > is the four fermion coupling constant. The chemical potentials of the two species can be written as 
/^i = /i + (5/i and /12 = ^ ~ (5/^, so that /i is the average of the two chemical potentials and their difference. The 
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Increasing Sfj. 




5|a, < 5|Xi 



5|a, > 5[Xi 



FIG. 3 (Color online). Pictorial description of the behavior of the Fermi spheres of two different populations of fermions, with 
up and down spin, in the presence of an attractive interaction and with increasing Sfi. For Sfi < Ao/V^ pairing takes place in 
the gray (red oline) region. For S/i > Aq/V^ the BCS homogenous pairing is no more energetically allowed and the system has 
a first order phase transition to the normal phase. 



effect of the attractive interaction between fermions is to induce the di-fermion condensate 

A(a;) 

{'4>six)'ipt{x)) = —^iia2)st , 



(9) 



which spontaneously breaks the global symmetry corresponding to particle number conservation. As a result the 
fcrmionic excitation spectrum consists of two Bogolyubov modes with dispersion laws 



e 



(10) 



with ^ = —fi+p^ /2m and Aq is the homogenous mean field solution. Without loss of generality we take 6fi > 0, then 



from Eq. ([TU]) we infer that tuning the chemical potential difference to values Sfi > Aq, the mode b becomes gapless. 
This phase corresponds to a superconductor with one gapped and one gapless fermionic mode and is named gapless 
homogeneous superfluid. 

In the above naive discussion we did not take into account that increasing Sfj, the difference between the free-energy 
of the superfluid phase, flsi and of the normal phase, f2„ decreases; eventually the normal phase becomes energetically 
favored for sufficiently large Sfi. In weak coupling it is possible to show tha t the two free-energies become equal at 
Sfii ~ Ao/\/2 (corresponding to the so-called Chandrasekhar-Clogston limit ( Chandrasekhari Il962 : Clogston . 1962I )). 
that is before the fermionic excitation spectrum becomes gapless. At this critical value of Sfi a first order transition 
to the normal phase takes place and the superfluid phase becomes metastable. The reason for this behavior can be 
qualitatively understood as follows. Pairing results in an energy gain of the order of Aq, however BCS pairing takes 
place between fermions with equal and opposite momenta. When a mismatch between the Fermi sphere is present it 
tends to disfavor the BCS pairing, because in order to have equal momenta, fermions must pay an energy cost of the 
order of Sfi. Therefore, when Sfi > cAq, where c is some number, pairing cannot take place. In the weak coupling 
limit, one flnds that c = l/\/2. This behavior is pictorially depicted in Fig. [3] 

Considering homogeneous phases, a metastable superconducting phase exists for Sfj, > Ao /y/2, but still the system 
cann o t develop fermion i c mass l ess modes, becatise at Sfi = An various ins t abilities appear ([Cubankova et al . 2010l 



20061: iMannarelli et all l2006at |Pao all l2006t ISheehv and Radzihovsk^ 120061: IWu and Yipl . l2003l) . I^Torder 



^ to 

explain what happens, let us consider the low energy spectrum of the system, which can be described considering 
the fluctuations of A(x) around the mean held solution Aq. The oscillations in the magnitude of the condensate 
are described by the Higgs mode, A(x), while the phase fluctuations are described by the Nambu-Goldstone (or 
Andcrsso n-Bogolyubov) mode (t>( x). Integrating out the fermionic degrees of freedom results in the Lagrangian 

m 



densities (Cubankova et al. 201 



C4, = A{dt(f>{x))' 



f(W(x))^ 



E . 



Cx = -CA(x)2 + D{dt(bix))^ - jiV^bix))^ 



(11) 
(12) 



The stability of the system is guaranteed when all the coefficients A, B, C, D, E are positive. A and D turn out 
to be always positive, then we define the three stability conditions 
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1. The Higgs has a positive squared mass: C > 

2. The space derivative of the Higgs must be positive: E > 

3. The space derivative of the NGB must be positive: B > 

For dfi > Aq the three conditions above are not simultaneously satisfied. Condition [1] is not satisfied because we 
are expanding around a maximum of the free-energy, and indeed a free-energy analysis shows that at — Aq the 
homogenous superfiuid phase becomes unstable. The fact that the condition [2] is not satisfied signals that the system 
is unstable toward space fiuctuations of the absolute value of the condensate, while condition |3] is not fulfilled when 
the system is unstable toward space fiuctuations of the phase of the condensate. Clearly the conditions [2] and |3] are 
related and tell us that when a large mismatch between the Fermi sphere is present, the system prefers to move to a 
non-homogenous phase. 

In the present toy-model the three conditions above are simultaneously violated in weak coupling at J/i/Ag = 1, 
but they are violated at different values of this ratio in the strong coupling regime. Moreover, with increasing coupling 
strength it is possible to force the system into a ho mogeneous gapless phase , but this happens when fi ^ — A, deep in 
the Bose-Einstein condensate (EEC) limit, see e.g. Gubankova et al. ( 2010t ). 



For 6fi < Aq, the homogenous BCS might also be energetically favored if there is a way of reducing Qg, and this is 
indeed what happens in the CFL phase where multiple interaction channels are available and the color and electrical 
neutrality conditions may disfavor the normal phase. This corresponds to the fact that C > and thus the condition 
[T] above is satisfied. However, the conditions [2] and [3] must be satisfied as well, but decreasing fig does not per se 
guarantee that those constraints are satisfied. Indeed, it is possible to show that in the weak coupling regime the 
gapless homogeneous phase is in general not accessible, because when (5^ > Aq, the solution with Aq 7^ is unstable 
with respect to the conditions [2] and [H 

Gauging the U{\) global symmetry, it is possible to show that the condition _B < is equivalent to the condition that 
the Meissner mass squared of the gauge field becomes negative, which corresponds to a magnetic instability. Therefore 
the magnetic instability is related to the fact that we are expanding the free-energy around a local maximum. This 
statement is rather general and indeed in Sec. III.FI we shall see that an analogous conclusion can be drawn for the 
2SC phase. Notice that increasing the temperature of the system does not help to recover from this instability 



(lAlford and Wangj . 120051) . Indeed, the effect of the temperature is to produce a smoothing of the dispersion law, 
which first has the effect of increasing the instability region to values Sfx < Aq. 

Summarizing, we have seen that for the simplest case of a weakly interacting two-flavor system, for S/j > Sfii the 
superfiuid homogenous phase is metastable, while for 6fi > Aq, fig it does not have a local minimum in Aq 7^ and 
it is unstable toward fluctuations of the condensate. In general, the three conditions above should be simultaneously 
satisfled in order to have a stable (or metastable) vacuum. The gapless phase is only accessible for homogenous 
superfluids deep in the strong coupling regime, for negative values of the chemical potential. 

A different possibility is that gapless modes arise at weak coupling in a non-homogenous superfluid. As we shall 
see in the following sections, the non-homogenous LOFF phase is energetically favored in a certain range of values of 
Sfi larger than the Chandrasekhar-Clogston limit, see Sec. Ill.Ci it is (chromo)magnetically stable, see Sec. III.FI and 
it has gapless fermionic excitations, see Sec. III.Dl It is important to remark that the p resence of a gapless fermionic 



It has gapless lermiomc excitations, see bee. lll.JJl it is important to remark tnat trie p resence or a gapless termiomc 
spectrum is not in contrast with the existence of superconductivity ( de Genned. 1966[l. e.g. type II superconduct ors 



have gapless fermionic excitations for sufficiently large magnetic fields (|de Genned " 19661 Saint- James et al. . 1969f ). 



B. Gapless 2SC phase of QCD 



The gapless 2SC phase ( g2SC in the following) of QCD was proposed in Shovkovv and Huand ( 20031 ) (see also 
Huang and Shovkovvl ( 20031 )) as a color superconducting phase which may sustain large Fer mi surface mismatches . 
Howev er, it was soon realized by the same authors that this phase is chromomagnetic unstable (jHuang and Shovkovv . 
2004ai rbh. meaning that the masses of some gauge flelds become imaginary. In the following we briefly discuss the 
properties of the g2SC phase at vanishing temperature, and then we deal with the problem of chromomagnetic 
instability. 

We consider neutral two-flavor quark matter at finite chemical potential. The system is described by the following 
Lagrangian density: 



£ = i'i-lf^d'^ ~m + /x7o) -0 + Ci„t , 



(13) 
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where ip = ipf , i = 1,2 a = 1,2,3 corresponds to a quark spinor of flavor i and color a. The current quark mass is 
denoted by m (we take the isospin symmetric limit m„ = = to), and Cint is an interaction Lagrangian that will 
be specified later. 

In Eq. (|13p . fi is the quark chemical potential matrix with color and flavor indices, given by 

2 

A* = fJ-ij.aP = (M Sij - [le Qij) Sap + --j=^Xs,{Ts,) apSij , (14) 

Qij = diag(QM, Qd) is the quark electric charge matrix and Tg is the color generator along the direction eight in the 
adjoint color space; ^e, ^J■a denote respectively the electron and the color chemical potential. Since fi is diagonal in 
color and flavor spaces, we can indicate its element with fiia, e.g. is the chemical potential of up blue quarks. A 
chemical potential along the third direction of color, ^3, can be introduced besides fig, but, for all the cases that we 
discuss in this section, we require the ground state to be invariant under the SU{2)c color subgroup; this makes the 
introduction of /X3 unnecessary. 

As interaction Lagrangian density we consider the NJL-like model 

Ant = Gs [(^AV')^ + (V'i75T-V')^] + Gd [{i^'^eeij5i')kj{i^eeij5i^^)k^] , (15) 

where ^p'-^ = Cip'^ denotes the charge-conjugate spinor, with C = 47270 the charge conjugation matrix. The matrices 
£ and e denote the antisymmetric tensors in flavor and color space, respectively; we used in the second term on the 
right hand side of Eq. (fT5|) the shorthand notation 

(V;'^eez75V')fc7 = [i^'ta^ijk^ap-tilbipjp) , (16) 

and an analogous expression for the other bilinear. In Eq. (1151) . two coupling constants are introduced in the 
sc alar-pseudoscalar quark-anti quark channel, denoted by Gs, and in the scalar diquark channel, denoted by Gd. 
In lHuang and Shovkovv (I2OO3I) the parameters of the model are chosen to reproduce the pion decay constant in the 



vacuum, /^r = 93 MeV, and the vacuum chiral condensate (uu)^/^ = {dd)^/"^ = —250 MeV. M oreover, an ultraviolet 



cutoff A is introduced to regularize the divergent momentum integrals. The parameter set of iHuang and Shovkovv 



(|2003[ ) is given by 

A = 653.3 MeV , Gs = 5.0163 GcV"^ . (17) 
The relative strength between the couplings in the quark-antiquark and quark-quark channels could be flxcd by 



a Fierz rearranging of the quark-antiquark interaction, see for example Buballal (l2005l ). The Fierz transforma- 



tion gives Gp/Gs = 0.75. However, non perturbativc in-medium effects might change this value. Therefore, in 
Huang and Shovkovv ( 2003') the ratio of Gd to Gs is considered as a free parameter. 



In lHuang and Shovkovv (200311 the authors consider only the case to = and vanishing chiral condensate. When 



TO 7^ the chiral condensate in the ground state does not vanish, but its effects are presumably negligible, giving 
a small shift of the quark Fermi momenta. This shift of Fermi momenta migh t change the numerical valu e of the 



electron chemical potential only of some few percent. Hence, the main results of iHuang and Shovkovvl (120031 ) should 
not change much if a finite value of the current quark mass is considered. 

Once the Lagrangian density is specified, the goal is to compute the thermodynamic potential. In the mean field 
(and one loop) approximation, this can be done easily using standard bosonization techniques. Neglecting the chiral 
condensate, the mean field Lagrangian density can be written within the Nambu-Gorkov formalism, see Sec. IIII.B.ll 
in the compact form 

/: = xt5-x-^, (18) 



where 



(19) 

is the Nambu-Gorkov spinor and the the gap parameter, A = A2sc£"^^ SijsG-f^, is included in the inverse propagator 



-1 _ f hf.d'' + /X70 



S-'^ A.,u - ^20) 

' AT ij^d'' - /X70 ' 
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as an off-diagonal term in the "Nambu-Gorkov space" . 

We shall focus on t he zero temperature r e gime (for a discussion of the rather uncommon temperature behavior 
of the g2SC phase see iHuang and Shovkovv which is relevant for astrophysical applications. The one loop 

expression of the thermodynamic potential can be determined from the inverse propagator in Ea. (j20[) : for vanishing 
temperature it is given by 



for a derivation see e.g. iBuballal (|2005l ). In the above equation fio is an irrelevant constant which is chosen in order 
to have a vanishing pressure in the vacuum. The second addendum corresponds to the electron free-energy (electron 
masses have been neglected). The last addendum is the contribution due to the quark determinant. The sum runs 
over the twelve fermion propagator poles, six of them corresponding to quarks and the other six corresponding to 
antiquarks quasiparticles: 

Ei,2 = \p\ T ^Jiub , (22) 
E3,4 = \p\ T fMb , (23) 



i?5,6 = d>+ V(|p|T/i)2 + A2gc , (24) 



Er^s = -Sfi+^Ji\p\Tf^? + Al^c, (25) 
and Eg — E^, Eio — Eq, Eh — Ei, E12 = Eg. Here we have introduced the shorthand notation 

M = M-Y + y, dn^—. (26) 
Using the explicit form of the dispersion laws, the previous equation can be written as 



127r2 127r2 127r2 27r2 



(\/(p-m)' + A2sc + \/(p + m)^ 



^2SC 



-26{5y. A2sc) p ^ (-5^ - ^{v--^?+^\^^ , (27) 



where m± = M i a/ Sy? — ^^^q. 



The value of A2SC is determined by the solution of the equation 



i9A2sc 



, (28) 



with the neutrality constraints. 



which fix the values of /ie a nd /xs. 



ns^-jr-^O, ng^-— =0, (29) 



The numerical analysis of iHuang and Shovkovv ( 2003 ) shows that Ms is much smaller than /ie and A2SC, both for 
A2SC > Sy and for A2SC < Si-i. As a consequence, it is possible to simplify the equations for the gap parameter and 
the electron chemical potential, ([28|) and ([29|) respectively, by putting fig = 0. Therefore , the properties of the s ystem 
depend only on the values A2SC and Me and on the couplings Gd and Gs. The result of lHuang and Shovkovv! (|2003l) 
can be summarized as follows: 

• For Gd/Gs ^ 0.8, strong coupling, the 2SC phase is the only stable phase 

• For 0.7 < Gd/Gs ^ 0.8, intermediate coupling, the g2SC phase is allowed for Sfi > A2SC 

• For Gd/Gs ^ 0.7, weak coupling, only unpaired quark matter is favored. 

In the g2SC phase the quasiparticle fermionic spectrum consists of four gapless modes and two gapped modes, 
whereas in the 2SC phase there are two gapless fermionic modes and four gapped fermionic modes. In the latter case 
the only gapless modes correspond to the up and down blue quarks that do not participate in pairing. 
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1. Meissner masses of gluons in the g2SC phase 

The diquark condensate of the 2SC phase induces the symmetry breaking pattern reported in Eq. ((5]); in particular 
the group SU{3)c CS" U{l)em is broken down to SU{2)c ® f/(l)em, where U{l)em is the gauge group corresponding to 
the rotated massless photon associated to the unbroken generator 



Q ^ Q cos ( 



-Tg sin^ 

e 



(30) 



where g and e denote the strong and the electromagnetic couplings respectively, Q is the generator of the electromag- 
netic gauge transfo rmations and Ts is th e eighth g enerator of the color transformations. The mixing coefficients have 
been determined in Afford et al. ( 2000bh (see also iGorbar ( 2000l )) and are given by 



cos 6 = 



sm f = 



The linear combination 



Tg = T8Cos6i+ -gsin6l, 
e 



(31) 



(32) 



is orthogonal to Q and gives the broken generator; the corresponding gauge field, which we shall refer to as the 8 
mode, acquires a Meissner mass. Actually, the NJL-like Lagrangian in Eqs.(fT3|) and (fT5|). has only global symmetries, 
but gauging the SU{2>)c group and the U{1) subgroup of SU{2)l x SU{2)r, one has that the spontaneous symmetry 
breaking leads to the generation of Meissner masses for the five gluons associate d to the broken g enerators. In order 
to compute these masses, we define the gauge boson polarization tensor, see e.g. Le Bellac ( 20001 ) . 



Tr[T>:^S{q)TlS{q-p)] 



where S(j)) is the quark propagator in momentum space, which can be obtained from Eq. (|20p. and 

r^^^5 7^ diag(T„-Tj) , 
r^ = e7^diag(Q,-Q) , 



(33) 



(34) 
(35) 



are the interaction vertex matrices. The trace in Eq. psp is taken over Dirac, Nambu-Gorkov, color and flavor indices; 
a,b ~ 1, . . . , 8 indicate the adjoint color and we use the convention that the component with a,b = 9 corresponds to 
the photon. 

The screening masses of the gauge bosons arc defined in terms of the eigenvalues of the polarization tensor and in 
the basis in which H^^ is diagonal the Debyc masses and the Meissner masses are respectively defined as 



M 



D.a 



iim n:^:i(o,p) , 

p— ^0 



A^L.a = -^lim(5.. + ^)n^i(0,p). 



(36) 
(37) 



Both masse s have been evaluated in the 2SC phase in iRischkd (j2000bf ): iRischke and Shovkovvl (j2002l ): and 



Schmitt et al 



(|2004l ). The Debye masses of all gluons are related to the chromo-electric screening and are al- 
ways real, therefore do not affect the stability of the 2SG and g2SC phases. The Meissner masses of the gluons with 

adjoint color a = 1, 2, 3 are always zero, because they are associated to the unbroke n color subgroup SU(2 )r.. 

For non- vanishing values of the Meissner screening masses have been evaluated in Huang and Shovkovvl (|2004al lbh. 
Gluons with adjoint color a = 4, 5, 6, 7, are degenerate and in the limit ^8 = their Meissner masses are given by 



MA 



^2SC 



2(5^2 



^2SC 



(38) 



The squared Meissner mass turn out to be negative not only in the gapless phase, A2sc/<5/i < 1, but also in the gapped 
phase, when A2sc/<^A* < V^- This result seems in contrast with the result of the previous section, where an imaginary 
Meissner mass was related to the existence of a local maximum of the free-energy arising at 5^ = A. However, from 
the analysis of the 2SC free-energy of the system, one can see that when A2Sc/<5/x = V2 the state with A2SC 7^ 
corresponds to a saddle point in the A — (5/i plane. The neutrality condition transforms this saddle point into a local 
minimum. However, as explained in the previous section, the gauge fields can be related to the fluctuations of the 
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gap parameter. These fluctuations can probe all the directions in the A — Sjj. plane around the stationary point and 
would result in a low energy Lagrangian with dispersion laws akin to the those discussed in Eqs. (|ll[) and (|12p with 
the coefficients E and B negative. 

Finally let us consider the 8 mode, associated to the broken generator defined in Eq. (|32l). The corresponding 
Meissner mass can be obtained diagonalizing the polarization tensor in Eq. ([33]) in the subspace a, 6 g {8, 9}, or more 
directly, by substituting Ta — )■ Tg '^^ the vertex factors of the polarization tensor. The squared Meissner mass of the 
8 mode turns out to be 

Ml,, = 'J^f^f^ ( 1 - -=^,(,, _ , (39) 



277r \ y/Sfi^ - a: 



2 

2SC 



and becomes negative for A2SC /SfJ- < 1. As shown in Gatto and Ruggieri ( 2007t) . the instability in this sector is 



transmitted to a gradient instability of the pseudo-Goldstone boson related to the U{1)a symmetry which is broken 
by the diquark condensate. Although in the vacuum the U{1)a symmetry is explicitly broken by instantons, at finite 
chemical potential instantons are Debye screened and U{1)a can be considered as an approximate symmetry, which 
is then spontaneously broken by the diquark condensate. 



C. The two-flavor crystalline color superconducting phase 

Since the homogeneous g2SC phase is chromomagnetically unstable, the question arises of the possible existence 
of a different superconducting phase for large mismatch betw een the Fermi sphere s. There are severa l cand idate 
phases, which include the gluonic phase (jGorbar et the solitonic phase (jNickel and Buballal . l2009t ) and 



the crystalline color superconducting phases. 

In this section we review some of the main results about the two-flavor crystalline phase. Firstly, we describe the 
one plane wave ansatz in the framework of the symple model discussed in Sec. III.Ai then, we turn to the crystalline 
color superconducting phase and report on the Ginzburg-Landau (GL) analysis of various crystalline structures. We 
relax the constraint of electrical and color neutrality, and treat the difference of chemical potentials between u and d 
quarks, 2Sfj,, as a free parameter. 



1. The one plane wave ansatz 



In Sec. lII.AI we have shown that in weak coupling the Chandrasekhar-Clogston limit (jChandrasekhaii Il962l : IClogston . 
1962h signals that the standard BCS phase becomes metastable, but does not forbid the existence of different forms of 
superconductivity. In particular, it does not forbid the existence of Cooper pairs with non- vanishin g total momentum . 



It was shown by Larkin and Ovchinnikov ([Larkin and Ovchinnikovl 11964) and Fulde and Ferrel (jFulde and Ferrell 



, in the context of electromagnetic superconductivity, that in a certain range of values of dfi it might be ener- 
getically favorable to have Cooper pairs with nonzero total momentum. For the simple two-level system discussed 
in Sec. III.AI Cooper pairs with momentum 2q can be described by considering the difermion condensate in Eq. (|9]) 
given by 



2iqx 



(40) 



and we shall call this state of matter the FF phase. 

Notice that this ansatz breaks rotational symmetry because there is a privileged direction corresponding to q. In 
the left panel of Fig. |4] the two Fermi spheres of fermions are pictorially shown and the red ribbons correspond to the 
regions in momentum space where pairing occurs. Pairing between fermions of different spin can only take place in 
a restricted region of momentum space and this implies that A < Aq. The reason why the FF phase is energetically 
favorable with respect to the normal phase, is that no energy cost proportional to dfj, has to be payed for allowing 
the formation of Cooper pairs. The only energetic price to pay is due to the kinetic energy associated to Cooper 
pairs. Actually, there is a spontaneous generation of a super-current in the d irection of q, which is balanced by a 
current of normal fermions in the opposite direction (jFulde and FerrellL [l96l . The gap parameter in Eq. (|40| can 
be determined solving a gap equation under the constraint that the modulus of the Cooper momentum, q, minimizes 
the free-energy. The result is that at S^i ~ (5/ii = Aq/v^ there is a first order phase transition from the homogeneous 
BCS phase to the FF phase. Increasing further 5^ results in a smooth decreasing of the gap function of the FF phase. 
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FIG. 4 (Color online). Pictorial description of the LOFF pairing in the weak coupling approximation. When 5/1 > Ao/v^ 
the BCS homogenous pairing is not energetically allowed, but pairing between fermions with total nonvanishing momentum 
can be realized. Left panel: In the FF phase pairing takes place in two ribbons on the top the Fermi spheres of up and down 
fermions, such that Pu + Pd = 2q, having opening angle 29q ~ 2 arccos((5^/g) ~ 67°, thickness A and angular width A/q, see 
Sec lII.E.T] Right panel: Structure obtained with two plane waves corresponding to two vectors 2q^ and 2q^ with relative angle 
(j>. The size and the opening angle of each ribbon is as in the FF phase. The structure with (j> ~ 180° is called the "strip". For 
illustrative purposes, we have greatly exaggerated the splitting between the Fermi surfaces, relative to the values used in the 
calculations reported in Sec. III. CI 



until at a critical value Sfi2 a second order phase transition to the normal phase takes place. In the weak coupling 
limit (5/i2 — 0.754Ao; the range [6^1,6^2] is called the LOFF window. In the LOFF window, the optimal value of q 
turn s out to be approxim ately constant, q « l.2dfi. 

In lAlford et al\ (I2OOII) the two-flavor QCD analog of the FF phase was presented. In this case the condensate 
has the same color, spin and flavor structure of the 2SC condensate, but with the plane wave space dependence 
characteristic of the FF phase, that is 



i2q-x aP3 

£ £ij3 



(41) 



As noticed in lAlford et 



(120011 ) ■ the FF condensate induces a spin-1 condensate as well; however, its effect is 
found to be numerically small, and it will be neglected here. As in the simple two-level system, in the two-flavor FF 
phase it is possible to determine the free-energy employing a NJL-like model with gap parameter in (|41l) . which can 
be determined solving the corresponding gap equation under the constraint that the value of q in Eq. (1411) minimizes 
the free-energy. The results are the same obtained in the two-level system; in particular the LOFF window and q 
have the same expressions reported above (but now 26 is the the difference of chemical potentials between u and d 
quarks). 



2. Ginzburg-Landau analysis 

From Fig.|4l it seems clear than an immediate generalization of the FF phase, can be obtained adding more ribbons 
on the top of the Fermi spheres, corresponding to different vectors q,™, with S {q}, where {q} is some set of vectors 
to be determined by minimizing the free-energy of the system and m is a label that identifies the vectors of the set. 
This in turn, corresponds to consider inhomogeneous CSC phases with a more general ansatz than in Eq. (j4ip . where 
the single plane wave is replaced by a superposition of plane waves, that is 

^^2^q■r _^ Aqn.e^*'?'" '- . (42) 

Assuming that the set of vectors {q} identifies the vertices of a crystalline structure it follows that at each set {q} 
corresponds a particular crystalline phase. We shall assume that the vectors q™ have equal length and then we can 
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write q"^ = qn'^, therefore the crystalhnc structure is determined by the set of unit vectors {n}. We shall also 
consider the simplified case that the coefficients A^™ do not depend on q"* and we shall indicate their common value 
with A. In other words, we shall consider condensates with 

p 

A(r) = A ^ g2i,«-.r- ^ ^43) 

rn—l 

where P is the number of vectors n™. The simplest example is clearly the FF condensate, depicted in the left panel 
of Fig.m characterized by a single plane wave, thus corresponding to P = 1. The case with P = 2 is reported in the 
right panel of FigUJ in this case the "crystalline" structure is completely determined by cfi, the relative angle between 
and n^; more complicated structures can be pictorially represented in a similar way. 

It is important to stress that the crystalline structure is determined by the modulation of the condensate, but the 
underlying fermions are not arranged in an ordered pattern, indeed fcrmions arc superconducting, that is they form 
a superfluid of charged carriers. 

The computation of the free-energy of a system with a general crystalline condensate cannot be obtained ana- 
lytically and also the numerical evaluation is quite involved. As a consequence, the use of some approximation is 
necessary. A viabl e method is the GL expansion of the free-energy, which is obtained expanding 51 in powers of A 



(jBowers and RajagopaL .20021) : 



2^./^.=P«A^ + |A^ + t^VC^(A«), (44) 

where the coefficients a, 13 and 7 are computed, in the one-loop approximation using, as microscopic model, a NJL 
model. The GL expansion is well suited for studying second order phase transitions but might give reasonable results 
for soft first order phase transitions as well. In the present case the expansion is under control for A/q ^ I and if 
the coefficient 7 is positive, meaning that the free-energy is bounded from below. 

For a given crystalline str ucture, the coefficients in Eg. (|44)) depend on Sfi and on the magnitude of q; the latter 



is fixed, in the calculation of Bowers and Raiagopall (2002), to the weak coupling value q = 1.26 ^i. For any value of 



Sfi the thermodynamic potential of a given structure is computed by minimization with respect to A and then the 
optimal crystalline structure is identified with that with the lowest free-energy. It is possible to compute analytically 
th e GL coefficients only for few s tructures: in general, they have to be computed numerically. We refer to the appendix 
of Bowers and Raiagopal ( 20021) f or details. 



In Bowers and Raiagopall (2002), twenty-three crystalline structures have been studied and among them, those with 
more than nine rings turn out to be energetically disfavored. This has been nicely explained in the weak coupling 
regime: in this case, as shown in the left panel of Fig. 2] for the FF phase, the pairing regions of the inhomogeneous 
superconductor can be approximated as rings on the top of the Fermi surfaces; one ring per wave vector n in the set 
{n}. The computation of the lowest or der GL coefficients shows that the intersection of two rings is energetically 



disfavored ([Bowers and Raiag opal 120021 ). As a consequence, it is natural to expect that in the most favored structure 



no intersecting rings appear. Since each ring has an opening angle of approximately 67 degrees, a maximum of nine 
rings can be accommodated on a spherical surface. 

This result can be quantitatively understood as follows. For the case of the two plane waves structure (right panel 
of Fig. |4]), in the weak coupling approximation there is one pairing ring for each of the two wave vectors. 

The quartic coefficient /3 depends on the angle </> between the two wave vectors and it diverges at 

(t)Q « 26'„ « 2 arccos — « 67° , (45) 

q 

which corresponds to the angle at which the two pairing rings are contiguous, meaning that for cf) < c/jq the two rings 
overlap. The latter case is energetically disfavored because, being /3 large and positive, the free-energy would be 
smaller. 

The divergence of the coefficient I3{(f>) at = (/jq is due to the two limits that have been taken to compute the 
free-energy, namely the Ginzburg-Landau and weak coupling limits. A detailed explanation of what happens will be 
given in Sec. IIII.BI when discussing a simple crystalline structure in the three-flavor case. In any case it is clear that 
the divergency of a GL coefficient means that the expansion is not under control, or more precisely, that the radius 
of convergence of the series (|33]) tends to zero. 

Among the structures with no intersecting rings, seven are good candidates to be the most favored structure. Within 
these seven structures, the octahedron, which corresponds to a crystal with P = 6 and whose wave vectors point into 
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the direction of a body-centered-cube (BCC), is the only one with effective potential bounded from below (that is, 
with 7 > 0). The remaining six structures, with P = 7,8 and 9, are characterized by a potential which is unbounded 
from be low, at least at the order A^. Even if in this case the free-energy cannot be computed, qualitative arguments 
given in lBowers and Raiagopal ( 2002 ) suggest that the face-centered-cube (FCC), with P = 8 and with wave vectors 
pointing towards the vertices of a face-centered-cube, is the fa vored structure. 

Of course, as the authors of Bowers and Raiagopal ( 2002f) admit, their study cannot be trusted quantitatively, 
because of the several limitations of the GL analysis. First of all, the GL expansion formally corresponds to an 
expansion in powers of A/q and therefore it is well suited for the study of secon d order phase transitions, but the 
condition that A/q 1 is not satisfied by all crystalline structures considered in iBowers and Raiagopall (|2002i ). In 
some cases the GL analysis predicts a strong first-order phase transition to the normal state, with a large value of 
the gap at the transition point. Moreover, it may happen that the local minimum for small values of A/q is not a 
global minimum of the system, as discussed in Sec. III.GI In this case the GL expansion in Eq. (|44p underestimates the 
free-energy of the system and is not able to reproduce the correct order of the phase transition. For a more reliable 
determination of the ground state one should consider terms of higher power in A, which are difficult to evaluate. 
Finally, the claimed favored crystalline structure, namely the face-centered-cube, has 7 < and a global minimum 
cannot be found unless the coefficient 0{A^) (or of higher order) is computed and found to be positive. 

Because of these reasons, the quantitative predictions of the GL analysis should be taken with a grain of salt. One 
should not trust the order of the phase transition obtained by the GL expansion and also the comparison among 
various crystalline structures may be partially incorrect, because it is not guaranteed that one is comparing the 
energies of the true ground states. 

On the other hand, the qualitative picture that we can draw from it, namely the existence of crystalline structures 
with lower free-energy than the single plane wave, is quite reasonable: crystalline structures benefit of more phase 
space available for pairing, thus lowering the free-energy. The symmetry argument is quite solid too, because it is 
based on the fact that configurations with overlapping pairing regions are disfavored, and as we shall see for one 
particular configuration in the three-flavor case in Sec. IIII.Bl one can prove that this statement is correct without 
relying on the GL expansion. Also, we shall show an interesting point, that the GL expansion underestimates the 
free-energy of the crystalline structures. And this happens not only in the presence of a global minimum different 
from the local minimum around which the GL expansion is performed, see Sec. III.Gi but also comparing the GL 
free-energy with the free-energy evaluated without the A/q expansion. 



D. Dispersion laws and specific heats 

The thermal coefficients (specific heat, thermal conductivity etc.) of quark matter at very low temperature are 
of fundamental importance for the transport properties and cooling mechanisms of compact stars. The largest 
contribution to thermal coefficients comes from the low energy degrees of freedom and it is therefore of the utmost 
importance to understand whether fermionic modes are gapped or gapless. Indeed, the absence of a gap in the spectrum 
of fermions implies that quasi-quarks can be excited even at low temperature and therefore the corresponding thermal 
coefficients are not suppressed by a factor « e^^l^ (which is distinctive of homogeneous BCS superconductors). 

In this section, we discuss the fermion and phonon dispersion laws in the two-flavor crystalline phases for low values 
of momenta. Then we employ the obtained dispersion laws for the computation of the specific heats. The results 
discussed here do not rely on the GL app roximation but are obtained by finding the zeros of the fu ll inverse propagator 
close to the nodes of the dispersion law (jCasalbuoni et al\ . 120031 : iLarkin and Ovchinnikov , 



1. Fermi quasi-particle dispersion law: general settings 

We consider a general di-fermion condensate A(r), and determine the quasi-particle dispersion laws looking at the 
zero modes of the inverse propagator of the system. Arranging the fields in the Nambu-Gorkov spinor as follows. 



the inverse propagator is given by 



^ -\ -e"/33£^,3A(r)* (5"^[Jy(£;-z-!;. V) + M^3).y] ' ' ^ ^ 
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where E is the quasi-particle energy and v is the Fermi velocity, that for massless quarks satisfies u = \v\ = 1. Then, 
the quasiparticle spectrum can obtained solving the eigenvalue equation 

{S-')rfx^^O- (48) 



Performing the unitary transformation 



= i^e'^'' "-"I" jC^, F°'= [e-'^" ""^ era) F°' , (49) 

it is possible to eliminate the dependence on (5/i in the eigenvalue problem and this corresponds to measure the energy 
of each flavor from its Fermi energy. The resulting equations for and Gf are independent of color and flavor indices, 
and therefore these indices will be omitted below. The eigenvalue problem reduces to solve the coupled differential 
equations: 

{E + iv ■ V)G - iA{r)F ^ , 

{E -iv ■V)F + iA{r)*G = . (50) 

These equations can be used to find the dispersion laws for any non- homogeneous A(r), and we shall consider 
here the periodic structures of the form given in Eq (j42p in order to determine whether there are gapless fermionic 
excitations. We shall prove that for any crystalline structure with real-valued periodic functions A(r) there exists a 
gapless mode iff the set {n} does not contain the vector n = 0. 

The proof is given below. Here we notice that this theorem does not apply to the case in which A(r) is not real, 
and indeed we shall show that the single plane wave structure does have a gapless mode. The theorem implies that 
any antipodal structure has a gapless mode, in particular the "strip" (corresponding to the structure depicted in the 
right panel of Fig. 2] for cj) = 180°) and the cube have gapless modes. On the other hand, the set of vectors which 
identify a triedral prism or a hexahedral prisms have a vector with n — 0, and the corresponding dispersion laws are 
gapped. 

Proof: 

For a periodic A(r), a general solution of the system in ([50]) is given by the Bloch functions 

G{r) = u(r)e''=-'' F{r) = w{r)e'''-'' , (51) 

where u(r) and w(r) are periodic functions. Notice that if k is real, the dispersion law is ungappcd, otherwise a gap 
is present in the excitation spectrum and the cigcnfunctions decrease exponentially with r. 

Taking E = m ([5D)l . the two solutions of the system of equations are given by F^ = G+ or by F_ = — G_, where 

G±(r) = exp(± A(r')^j , (52) 

with ry = r • v. 

Comparing this expression with Ea.([5T|). it is clear that theses solutions corresponds to Bloch functions with 
Re(/c) = 0. If A(r) has a term with n = 0, it means that A(r) = A + /(r), where A is a constant and /(r) is a 
periodic non-constant function. In this case Eq. (|52p has an exponential behavior of the type G(r) oc exp (±Ar) and 
therefore the spectrum is gapped. On the other hand, if in the expansion of A(r) no term with n = is present then 
the imaginary part of k vanishes and the spectrum is gapless c.v.d. 

The fcrmion dispersion law for the gapless modes can be determined using degenerate perturbation theory for 
^ = ]3 — /i <C 9, i.e. close to the Fermi sphere. At the lowest order in S^/q one finds that 

E{v, = , , = c{vK , (53) 

where c{v) is the velocity of the excitations, and 



A±{v) = — f drexp ±2 / A(r') 



(54) 



where Vc is the volume of a unit cell of the lattice. The energy of the fermionic excitations depends linearly on the 
residual momentum ^, but the velocity of the excitations is not isotropic. 
Let us now specialize ((53)) to the case of the strip, 

A(r) = 2Acos(2q • r) , (55) 
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FIG. 5 (Color online). Velocity of the fermionic quasiparticles as a function of 2 = cosi? for 2/S./{qv) = 0.1 (solid red line), 
2A/{qv) = 0.2 (dashed green line), 2A/{qv) = 0.3 (dotted blue line). 



which corresponds to the condensate in Eq. (|43)) with m ~ 2 and n? — —n} — n and strictly speaking docs not 
describe a crystal, but a condensate that is modulated in the n— direction. The coefficients in the dispersion law are 
given by 

^(strip) _ ^(.trip) ^ j^i^^ ^ (5g) 

where Iq{z) is the modified Bessel function of the zeroth order. Therefore the velocity of the fermionic quasiparticles 
has the analytic expression 

Cstrip(t') = _ /^^^ , (57) 



q-v 



which has the important property to vanish when v is orthogonal to n. The reason is that, in the direction orthogonal 
to n the gap is constant and its effect is equivalent to a potential barrier. Taking n = (0, 0, 1), the dispersion law is 
symmetric with respect to rotations around the z— axis, for inversions with respect to the plane z = and depends only 
on the polar angle between v and the z— axis. In Fig. [5]we report a plot of the velocity of fermionic quasiparticles 
as a function of cosd, for three different values of the ratio A/{qv). In the ultrarelativistic case, w = 1, the relevant 
case is A/q < 1 and we see that the dispersion law of fermionic quasiparticles is not much affected by the condensate 
for cos > 0.2A/q' and it is the same of relativistic fcrmions. On the other hand for small values of cos the fermionic 
velocity is exponentially suppressed and vanishes for A = 0, meaning that fermionic quasiparticles cannot propagate 
in the x — y plane as discussed above. 

The fact that the dispersion law is linear in ^ for small values of the momentum does not assure that it is linear 
for any value of the momentum. Considering <^ 1, it is possible to solve th e Eq. (|50|) for k in ([5T|) along the 



z— direction, without restricting to low momenta (jLarkin and Ovchinnikovl . 119641 ): the result is that 



meaning that the dispersion law is linear in the residual momentum, only for k/q <^ 1. 

For the octahedron, whose six wave vectors point into the direction of a BCC structure, the corresponding gap 
parameter can be written as 

A(r) = 2A[cos(2(7a:) + cos(2gy) + cos{2qz)] . (59) 
It is easy to show that in this case the integral in Eq. ([M)) factorizes, and the dispersion law is gapless with velocity 

CBCc(^') = — r~r\ 7~T\ / . \ • (60) 
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FIG. 6 (Color online). Velocity of the fermionic quasiparticles in the BCC crystalline structure (left panel) and in the FCC 
crystalline structure (right panel) as a function of the polar angles i? and (f, for 2A/{qv) = 0.3. 



The corresponding plot is reported in FiglHl left panel, where the unit vector v has been expressed by the polar angles 
1? and Lp. The plot has been obtained for 2A/{qv) = 0.3 and considering -d G [0, tt] and ip G [0, 2tt]. Note that according 
with Eq. (j60p . the velocity of the fermionic quasiparticles vanishes along the planes Vx = 0, Vy = and Vz = 0. 

For more complicated crystalline structures it s not possible to have an analytic expression of the fermionic velocity, 
one notable example is the FCC structure which is defined by the condensate 

4 

A(r) = 2A ^ cos(2g n" • r) , (61) 



where = y i(l, 1, 1), ~ y i(l, 1, — 1), n'^ = y i(l, — 1, 1) and = y ^(—1, 1, 1). Upon plugging this expression 
in Eq. ([M)) we obtain that 

^iFCO^fgyf rfyexpli^sj, (62) 



-l"^ J cell I 

where the integration is over the elementary cell of volume {tt /qf" and 

B = VSw { ^'^^^ + y + z) ^ sin 2q{x + y - z) ^ sin 2q{x - y + z) ^ sin 2q{-x + y^^) \ ^^^^ 



and it is then easy to show that A^^^^^ = A^^^^'^ = Af^cc), 

Expressing the components of the unit vector t> in _B as functions of the polar angles § and Lp, and upon substituting 
in one has that the quasiparticle velocity 

CFCc(^,y>) = ^(PCC)(^^^) ' (64) 

has the behavior reported in Figl6l right panel. The plot of cfcc(i?j<p) has been obtained for 2A/{qv) = 0.3 and 
considering d £ [0, tt] and ip £ [0, 27r]. The velocity of the fermionic quasiparticles is equal to 1 almost everywhere with 
the exclusion of a restricted region which corresponds to the zeros of B{d,tp), given by the solutions of v ■ n"^ = 0. 

In conclusion, we have shown that for various crystalline structures the fermionic spectrum is gapless. For k/q <ti 1 
the dispersion law is linear in momentum and the fermionic velocity vanishes along the planes orthogonal to the 
direction of the vertices of the reciprocal lattice. Th is result remains valid als o for massive quarks, since the effect of 



the quark mass can be accounted for by reducing v (jCasalbuoni et all l2002a[ ) 
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Regarding the FF condensate we take the direction of the Cooper pair total momentum 2q along the z— axis. In 
this case, the quasi-particle spectrum has the analytical expression 

E± = qv, ± + A2 , (65) 

where the quasi-particle energies are computed from the corresponding Fermi energies Hu.d- Eq. ()65|) is the dispersion 
law of quasi-particle {E± > 0) or hole states {E± < 0). As for the case of the strip and of the FCC, the dispersion 
law depends on the direction in coordinate space that one is considering. However, contrary to what happens for real 
crystalline structures, in the FF phase there are directions along which the dispersion laws are gapped and directions 
along which the dispersions law is gapless. 

The zeros of the dispersion laws ares located, for a given value of cos??, at 

ICol = \/((?i;cost?)2- A2 , (66) 

and from Eq. (p5)) one can see that E_ is gapped for cost? > A/qv, while E^ is gapped for cosi? < —A/qv. Therefore 
for (gwcosi?)^ < A^, both dispersion laws would be gapped but in this case the zeros of the dispersion law becomes 
imaginary, meaning that both dispersion law cannot be simultaneously gapless. Note that for A/qv > 1 the gapless 
modes disappear. 



2. Specific heat of the Fermi quasi-particles 

The contribution of the Fermi quasi-particles to the specific heat per unit volume is given by 

j 

n(Ej, T) is the Fermi distribution function and the sum is over all the fermionic modes. Considering the low temper- 
ature range, T << A, which is relevant for astrophysical applications, and using Eq. (j53p one obtains for a generic 
two-flavor crystalline structure 

where the second t erm on the r.h.s is the c ontribution of unpaired quarks. This expression can be evaluated in closed 



form for the strip (jCasalbuoni et al\ . 120031 ) 



(strip) _ V^J' . P. / 1 /O. 1 /O 1 . r A /„„^2^ ^ 



1F2 (-1/2; 1/2, 1; {A/qvf) + ^ , (69) 



where 1F2 denotes the generalized hyperg eometric function ( Gradshtevn and Rvzhikl . Il980l ). Differently from the 



-nalysis of iLarkin and Ovchinniko here v is not small and we can take A/qv — ?► near the second order 



ar ^ 

phase transition. Since for small A/qv one has iF2(— 1/2; 1/2, 1; [A/qvY) ~ 1 — {A/qv)'^, it is easily seen that the 
normal Fermi liquid result is obtained for A = 0. On the other hand, for non-vanishing A, the specific heat turns out 
to be smaller. 

In the case of the FF state, the dispersion law of the quasi-particles is given by ([65|) and using Eq. (|67| one obtains 
that in the small temperature limit (T ^ A) and for A < q 



(FF) A2 2^l\ 

In the above equation, the first and the second addenda correspond to the contribution of the paired and the unpaired 
quarks respectively. The paired quarks contribution depends linearly on temperature because the quasi-particle 
dispersion law (p5|) gives rise to gapless modes when A/qv < 1. There is also a contribution to the specific heat that 
comes from gapped modes, but this contribution is exponentially suppressed with the temperature. 
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3. Effective Lagrangian of phonons and contribution to the specific heat 



The low energy spectrum of a periodically modulated condensate, besides Fermi quasi-particles consists also of 
masslcss NGBs which originate from the spontaneous b reaking of rotational sym metries. Note that for two-flavor 



QCD in any crystalline phase, these are the only NGBs ( Casalbuoni et al . 2002bl lc). because no global symmetry of 



the system is spontaneously broken. Since the modulation of the condensate is associated to a crystalline structure, 
the NGBs describe the vibrations of the crystal and are for this reason called phonons. These phonons are not 
associat ed to pressure oscilla. tions, rather it can be shown that they are related to chemical potential oscillations. 



see e.g. 



Anglani et al. 



(j201l[ ). For any crystalline structure one can deduce the general expression of the phonon 
Lagrangian from the symmetries of the system, but the various coefficients have to be evaluated by a microscopic 
theory. These coefficients are related to the elastic properties of the crystalline structure and we shall present an 
analysis of the shear modulus of various three-flavor crystalline structures in Sec. IIII.DI 

In the FF phase there is one single phonon field, 4>, and one privilege d direction corres p ondin g to q. The leading 



order (LO) Lagrangian density in the momentum expansion obtained in lCasalbuoni et all (|2003l ) is given by 



1 



(71) 



where V|| = n • V, Vj^ = V — nVn and n is the unit vector parallel to q. The breaking of the rotational symmetry 
in the underlying microscopic theory implies that the velocity of propagation in the direction parallel to n, v\\ is 
di fferent from the v elocity of propagation in the direction orthogonal to n, v±. These velocities have been computed 
in 



Casalbuoni et al, (2002c,l employing as microscopic theory a NJL-like model and the result is 



vu = cos^ 0, 



v^, = - sin^ ( 
2 



with cosOq ~ Sfi/q. The dispersion law, relating the phonon quasi-momentum k and energy tu, therefore is 



The contribution of phonons to the specific heat at small temperatures is given by 

S^2 



JFF) 
-V 



15w i v\ 



■ T (phonons) 



(72) 



(73) 



(74) 



Since the strip has the same space symmetries of the FF phase and one phonon field, the Lagrangian density has 
the same expr ession reported in Eq.(|71|). but the values of the longitudinal and transverse velocities are different 
(|Casalbuoni et al. 20023). Thus, for the strip the expression of the specific heat is formally the same reported above. 

The low energy oscillations the FCC crystal structure are described by three different modes and the LO Lagrangian 
density for the three phonon fields is given by 



1=1,2,3 



1=1,2,3 



4=1,2,3 



i<j"=l,2,3 



(75) 



where a, h and c arc three coefficients to be determined by the microscopic theory. Notice that the Lagrangian is 
rotationally invariant, therefore it is more symmetric than the underlying theory. The reason is that phonons are 
long wavelength fluctuations and therefore for the cubic structure they are not sensitive to the modulation of the 
condensate. This is different from the FF or strip structure, because in that case the microscopic symmetry is the 

same symmetry of the macroscopic level. 

The microscopical calculation of the coefficients in ([75|) has been done in ICasalbuoni et employing a 

NJL model, and the result is 



1/12 



, 



(3 cos^ 



1)/12. 



(76) 



The corresponding spcciflc heat contribution has been evaluated numerically in ICasalbuoni et al\ (j2003l ) , and the 
result is 



„(FCC) 



(phonons) . 



(77) 
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E. Smearing procedure 



In Section [11.0.21 we have discussed the Ginzburg-Landau expansion of the free-energy for various CCSC phases. 
Since the GL expansion has several hmitations and is rehable in a hmited number of cases, 

it is useful to derive a different approximation scheme. In the previous section we have presented an approximation 
that allows to deal with the dispersion law of fermionic quasip articles close to th e gaple ss momentum. Here we shall 
discuss the smearing procedure, an approximation developed in Casalbuoni et al. ( 20041 ) which is valid when A/Sf^i is 
large. 



1. Gap equations 

In order to define the smearing procedure wc first consider the FF phase. Although this case can be solved exactly, 
it is useful to consider it to fix the notation and introduce some definitions to be used later on. We shall consider 
Cooper pairing of the massless up and down quarks, chemical potential fiu, fJ-d, and we define fj, = {fiu + fJ'd)/'2 and 
SfJ- = \fJ-u - IJ-d\/ 2 < /X. 

As detailed in ICasalbuoni and NarduUil (|2004l ). fr om the mean field Lagrangian of the FF state 

Caond = -^e"'''e.j3(CV'f A(r-) + c.c.) + (L ^ R) - -A{r)A*ir) , (78) 
one can deduce the zero temperature gap equation, which in the HDET approximation is given by 



An 



where p ~ is the density of states in two-fiavor QCD and 5 is the ultraviolet cutoff taken equal to /i/2; moreover 
it is assumed that 6 ^ q ^ Sfi. The quasi-particlc dispersion laws have been obtained in Eq. (|65p . but energies are 
now measured from the common energy level ^, and therefore 



Eu,d = ±SfiTq-v + + A2 . (80) 

In Eq. (|79p the contribution of hole excitations is taken into account simply multiplying the contribution of quasi- 
particles times two. The reason is that the difference between the contribution of quasi-particles and quasi-holes is of 
order Sfx/fi. 

We observe that in general 

1 - 9(-x) - e{-y) = 9{x)0{y) - d{-x)e{-y) , (81) 
and since and Ed cannot be simultaneously negative, then 

1 - e{~Eu) - ei-Ed) = 0iEu)0iEd) , (82) 
so the integration in Eq. (|79p is over a restricted region named the pairing region (PR), defined by 

FR= {{^,v-n)\E.u > OandEd > 0;^ < 5} . (83) 
More explicitly, the pairing region is defined by the condition 



Max <^ -1, Zg - } < vn < Min n , + \ , (84) 



with 



= ^ . (85) 

q 

From the above definition it follows that for small values of A the pairing region is centered at i? = arccoszg, has an 
angular width of order A/q and a thickness of order A, corresponding to the upper red ribbon in the left panel of 
Fig. m The lower red ribbon can instead be obtained considering hole excitations and the same reasoning used above. 
In the limit (S/i/^ — the two ribbons have exactly the same dimensions. 
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Thus. Eq. (|79|) can be written in a different way: 



gp f f dQ, A gp f dn A^s 



2 J JpR 4^ Ve+^ 2 7 47r 7o VF+A^ 



A=fy y^^ ^rfc^=== = ^ / ^ / ..r".., ^ (86) 

where Aoff = Aoff(v • n,^) is defined as 



A for v-n) ePR 
A^,=AeiE..)0iE,) = i^^ elsewhere. 

The above procedure defines the smearing procedure for the FF phase; it can be extended to the case of P plane 
waves, Eq. (|43p . generalizing the results of the previous equations, assuming that in the mean field Lagrangian one 
can substitute A with Ae {v,£o), where 

p 

AB(t;,4) = 5lAcff(^-"'"'^o) , (88) 

n=l 

meaning that A^; = 71 A, where n = (1, . . . , P). We can thus generalize the pairing region to 

Vn = {{v,0\^E{v,e)^nA}. (89) 

Note that in this equation we have made explicit the dependence on the energy io instead of that on the residual 
momentum, because the pole position is in general in 



= + ^'A2 . (90) 



Correspondingly, the gap equation is 



PA^i^ — / (91) 



.gp r dn r dipdj Asivjo) 

2 J An J 27T £l-e-Al(v,e„) ' 

which generalizes Eq. (|86|). The origin of the factor P on the l.h.s. of this equation is as follows. The Lagrangian 
contains the term 

A^rOAW ^^^^ 
9 

which, when averaged over the cell, gets non vanishing contribution only from the diagonal terms in the double sum 
over the plane waves and each plane wave gives a separate contribution. 

The energy integration is performed by the residue theorem and the phase space is divided into different regions 
according to the pole positions. Therefore the gap equation turns out to be given by 

, 2(5 ^ f f dn ,^ nA 



where we have made use of the equation 



2 2^ 

-=ln— , 94 

gp Ao 



relating the BCS gap Aq to the four fermion coupling g and the density of states. The first term in the sum, 
corresponding to the region Vi, has P equal contributions with a dispersion rule equal to the Fulde and Ferrel case. 
This can be interpreted as a contribution from P non interacting plane waves. In the other regions the different plane 
waves have an overlap. Since the definition of the regions Vn depends on the value of A, their determination is part 
of the problem of solving the gap equation. 

Said in a different way, in the smearing procedure the dispersion relation of the quasi-particles has several branches 
corresponding to the values nA, n = 1, • • • , P. Therefore, the following interpretation of the gap equation (|93p can 
be given. Each term in the sum corresponds to one branch of the dispersion law, i.e. to the propagation of a gapped 
quasi-particlc with gap nA, which is defined in the region Vn- However, the regions Vn do not represent a partition 
of the phase space since it is possible to have at the same point quasi-particles with different gaps. 
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FIG. 7 (Color online). Free energies of the octahedron (dashed red line) and of the FCC (solid blue line) crystalline structures as 
a function of (J^/Aq. The octahedron is the favored structure up to Sjj, « 0.95Ao; for 0.95Ao < Sfi < 1.32Ao t he FCC is favored. 
Here, for each value of Sfi, the values of Zq and A are those that minimize the free-energy. Adapted from ICasalbuoni et all 



2. Numerical results: free-energy computation 



In Table U we report the results obtained at the Chandrasekhar-Clogston point ((5/.( = Sfii) for four crystalline 
structures, respectively the FF (P = 1), the strip (P = 2), the octahedron (BCC) (P = 6), and the FCC (P = 8). 
The table shows that, among the four considered structures, the favored one at = dfii is the octahedron, which 
however does not have the largest gap A. The gap parameter determines the extension of the pairing regions, see e.g. 
Eq. (|84p . but a free-energy gain may result from having many small non-overlapping pairing regions, as well. Indeed, 
the strip has the largest gap parameter, but not the largest free-energy, presumably because the pairing occurs only 
in four large ribbons (two for each vectors q'", sec the right panel of Fig|3]). Indeed, as we have shown in the previous 
section, the width and the thickness of the ribbons arc proportional to A. On the other hand, for the octahedron 
the pairing occurs in smaller ribbons, but there are twelve pairing regions, which give a large contribution to the 
free-energy. 

For the case (5/i ^ (5/ii, we report in Fig. [7] the plot of the free-energy of the octahedron (dashed line) and of the 
FCC (full line) structures as a function of S^/Aq. The octahedron is the favored structure up to (5/i w 0.95 Ag; for 
larger values of the chemical potential mismatch the FCC structure is favored, up to 6fj, ~ 1.32Ao where the normal 
phase becomes favored. Also in this case we find that the FCC gap parameter is smaller than the octahedron gap 
parameter, but in the FCC there are more pairing regions than in the octahedron. 

In Table |TT] we report numerical results for S112 for each crystalline structure and the computed order of the phase 
transition between the crystalline phase and the normal phase. We have reported also the values of Zq, see Ea.(|85l). 
and of the discontinuity in A/Aq at Sfi = Sfi2 — 0^. 

In the smearing approximation both the order of the tr ansition and the point where the transition occurs are 
different from those obtained within the GL approximation (jLarkin and Ovchinnikovll 19641 ). However, the difference 
in Sfi2 is ^ 10% and in Zq is 17% (in GL for any structure Zq ss 0.83). On the other hand, in agreement with 



p 


Zg 


A 
Ao 


2(1 


1 


0.78 


0.24 


-1.8 X 10"^ 


2 


1.0 


0.75 


-0.08 


6 


0.9 


0.28 


-0.11 


8 


0.9 


0.21 


-0.09 



TABLE I The gap, Z g — S/i/q and the free-energy at Sfi — Sfii — Ao/\/2 for different crystalline structures. From 
ICasalbuoni et a/.i l|2004l ) . 
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the GL results, the structure with 6 plane waves is more favorable than the structure with 2 plane waves or with 
1 pane wave. Therefore, increasing the pairing region leads to an increase of the free-energy (in absolute value) 
making the configuration more stable. An interesting result of the smearing procedure is that with increasing 6fj, the 
configuration with 6 plane waves is superseded by the structure with 8 plane waves, which has a lower A but a larger 
number of pairing regions. This has to be contrasted with the result of the GL approximation, see Sec lII.C?2l where 
the octahedron is always disfavored with respect to the FGC (but one should consider that the GL free-energy of the 
FCC is not bounded from below). 



F. Chromomagnetic stability of the two-flavor crystalline phase 



As discussed in Section llLB.ll the 2SC phase is chromo- magnetically unstable for Sjj > Sfii = Ao/V2- Further in- 
creasing (5/x one eventually enters into the g2SC phase which is chromo-magnetically unstable as we ll. This could be in- 
terpret ed as the tendency of the system to generate a net momentum of the quark pair, as shown in lGiannakis and Ren 
(l2005al) . Thus the chromo-magnetic instability can be interpreted as as a tendency to develop quark current, which 
in turn is equivalent to the FF phase, where diquark carry momentum 2q. 

In this section, we review the results of iGiannakis and RenI (|2005a[) . which relate the Meissner mass of the 8 mode 
of the 2SC to the momentum susceptibility and the n the computation of the M eissner tensor in the crystalline phases 
( Giannakis et ai . 2005 : Giannakis and Ren . 2005al lbl: Reddv and Rupak . 2005 ). 



1. Momentum susceptibility 

The response of the thermodynamic potential of the 2SC phase to a small momentum, 2q, of the quark pair can 
be computed absorbing the phase of the condensate into the phase of the quark fields; the net effect is a shift of each 
quark momentum by q, see Sec. IIILBI for more details on this procedure. Therefore, the one loop effective action (in 
presence of background gauge fields) can be computed using the same steps that lead to Eq. (pij) . Expanding the 
thermodynamic potential around q = one has, at the lowest order, 

n^n2sc + l)Cq\ (95) 

where the momentum susceptibility is given by 

^ = ^ E / Tr [r.5(p)r,5(p)] , (96) 

where S(jj) is the quark propagator in momentum space and T = diag(7, —7) is the appropriate vertex factor. 
Notice that we are considering an expansion of the thermodynamic potential for small q, therefore the momentum 
susceptibility does not depend on q, meaning that S{p) is the 2SC quark propagator which can be obtained inverting 
the expression in Eq. (|20p . The momentum susceptibility has an expression similar to the space component of the 2SC 
polarization tensor, see Eq p3|) . at vanishing momentum. In particular it is possible to show that it is proportional 
to the squared Meissner mass of the 8 mode. The reason is that the only non-vanishing contribution to both the 
momentum susceptibility and to the Meissner mass of the 8 mode are determined from the red-green color sector 
(blue quarks do not carry a condensate and do not mix with red and green quarks). But in this color sector Tg is 
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0.754 


II 


0.83 
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0.83 


I 


1.0 


0.81 
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L22 
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0.95 


0.43 
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L32 


I 


0.90 


0.35 



TABLE II The values of 5fi2, Zq = Sjj,/q, the discontinuity of A/Ao and the ord er of the phas e transition between the CCSC 
phase and the normal phases for different crystaUine structures. From lCasalbuoni et al\ l|2004f ). 
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proportional to the identity and thu s the co rresponding vertex factor is proportional to T. A more detailed discussion 
can be found in iGiannakis and RenI (l2005a|), where it is shown that 



(97) 



where the mass of the 8 mode in the 2SC phase is given in Eq. (|39|) . 

At the transition point between 2SC phase and the g2SC phase, = A, one finds Af|^ g ^^"^ thus, Eq. ([57)1 
implies that the system is unstable towards the formation of pairs with nonvanishing net momentum; indeed a negative 
JC in Eq. (|95|) implies that there is a gain in free-energy if q 7^ 0. Therefore, the chromomagnetic instability of the 8 
mode leads naturally to the FF state. 

Before turning to the computation of the Mcissner tensor in the crystalline phase, we comment briefly on the 
absence of total currents induced by the net momentum of the quark pair. The value of the total momentum is 
determined minimizing the free-energy, thus dQ/dq ~ 0, and the stationarity condition is equivalent to the vanishing 
of the tadpole diagram. 



(?/'7'0) = , 



(98) 



which implies that no ba ryon matter current is generated in the ground state. Analogously, one can show 
( Giannakis and Ren . 2005al) that electric and color currents vanish as well. In particular, the residual SU{2)c'SiU{1)q 
symmetry implies that {-tpTAilj) = for A 8, with Ta defined in Eqs. ([34|) and ([35]) : moreover, for A = 8 one finds 



Js = 



(99) 



which vanishes because of the minimum condition, Eq. (|98p . Therefore, no total current is generated in the FF state. 



2. Meissner masses in the FF phase 



Since the FF phase has the same gauge symmetry breaking pattern of the 2SC phase, it has five massive gluons . 
The computation of the Meissne r tensor of gluons in the two- flavor FF phase has been done in lCiannakis et al\ ( 2005 ) 
and IGiannakis and Ren (l2005bl) . neglecting neutrality conditions and considering the isospin chemical potential, 
Sfi = ^e/2, as a free parameter. The Meissner masses of the gluon fields with adjoint color a = A, ... ,7 are all equal 
and can be written as 



AS, 



r,2 



B 



(100) 



where we have decomposed the Meissner tensor into longitudinal and transverse components with respect to q. The 
Meissner tensor of the 8 mode can be decomposed in a similar way 



(101) 



The coefficients A and C are called the transverse Meissner masses; similarly, B and D are the longitudinal Meissner 
masses. As discussed in the previous section, the mass of the rotat ed eighth gluon is relat ed to the variation of the 
free-energy with respect to q. In the FF phase it has been shown by IGiannakis et al. ( 2005 ) that the precise relation 
is the following 



C = — 
12 



D ^ — 
12 



.9 +^ 



.9 +^ 



1 dn 

W\W\ 

3 J ' 



(102) 
(103) 



If the phase with q 7^ is a minimum of the free-energy, then both the conditions dil/d\q\ = and d'^rt/d\q\'^ > 
must be satisfied. Thus, C = and D > at the minimum. As a consequence, the Meissner tensor of the 8 mode is 
purely longitudinal and positively defined. 
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Th e coefficients A, B a nd D can be computed analytically. We refer the interested rea der to the origi n al ar- 
ticle ( Giannakis and Ren . J2005bl) . Here, it is enough to consider the small gap expansion (jCiminale et ali l2006t 
Giannakis and Ren . l2005bl) : in this approximat on scheme one has 



A4 



B = 



D = 



967r2 5^4(^-2 
A2 



1)2 



8^2 5^1^{Zq'-l) 



A2 



3.9V <5A*'(2g-' - 1) 



(104) 
(105) 
(106) 



The message of the above equations is that the Meissner tens or is positively defined for t he one plane wave ansatz, 
within the small gap parameter approximation. The authors I Giannakis and RenI (|2005br ) argue that for a multiple 
plane wave structure, the situation will be better (within the small gap expansion). As a matter of fact, at the jS^i} 
order, the Meissner tensor is purely longitudinal, and the longitudinal components arc positive at the minimum of 
the free energy. Therefore, if the expansion in plane waves contains at least three linearly independent momenta, the 
Meissner tensor will be positive definite, being it addi tive with respect to d ifferent terms of the plane wave expansion 
to order A^/ ^/x^. This has been explicitly checked in Ciminale et al. ( 20061 ). 

Besides, in iGiannakis and Ren ( 2005b ) the numerical computation of the Meissner masses, beyond the small gap 
expansion, has been performed. The results can be summarized as follows. First of all, the FF phase is found to be 
more stable than the 2SC phase in the range 



0.706 Ao <6ii< 0.754 Ao . 
Within this range, the LOFF gap is within the range 

< A < 0.242 Ao , 

and the longitudinal Meissner masses, B and D, turn out to be positive within the range 

< A < 0.84 Ao ; 

moreover, the transverse mass A of the gluons with a = 4, . . . , 7 turns out to be positive within the range 

0<A<0.38Ao. (110) 



(107) 
(108) 
(109) 



Since the LOFF window for the FF state (|108|) is contained in both the intervals (|109p and ()110p . the FF state is free 
from the chromomagnctic instability, as long as it is energetically favored with respect to the homogeneous phase. 



G. Solitonic ground state 



The analysis of various crystalline phases h as shown that a periodic structure is energetically favored for mismatched 
Fermi spheres. In iNickel and Buballal (|2009f ) a generalization of the crystalline structure has been proposed in order 
to explore whether the more complicated periodic condensate (we do not report the color-flavor structure because it 
is the same of Eq. (|4T|)): 



A(z) = J2 e' 



2ikqz 



(111) 



may be energetically favored with respect to standard crystalline structures. The wave vector q is taken along the 
z-axis, thus the condensate corr esponds to a band s t ructu re along the z— direction. Assuming that the condensate is 
real (which is the assumption of Nickel and Buballal ( 20091 )). one has Ag ^ = A* and Eq. (|111|) can be rewritten as 



A(z) = 2 ^ A„ cos(2ngz) 



(112) 



Written in this form, it is clear that the ansatz is a generalization of the strip structure and amounts to consider 
higher harmonics contributions to the gap function, which are not included in the strip. 
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0.72 0.74 



FIG. 8 (Color online). Left panel: profile for the gap function, at several values of q, obtained from the self-c onsistent solution of 
the ga p equation, at S/j, = O.TAbcs- Right panel; physical wave vector magnitude as a function of S/j,. From lNickel and Buballal 
(|2009h . 



For any value of Sfi, the ground state in Nickel and Buballa ( 2Q09[) is determined using a two-steps procedure. 
Firstly, the magnitude of q in Eq. (|112p is fixed, and the profile A(z) is determined by solving the gap equation. 
Then, among several choices for q, the physical value corre s pond s to the one which minimizes the free-energy. In the 
left panel of Fig. El which is taken from Nickel and Buballal ( 2009f ). we plot the profile A{z) obtained by the numerical 
solution of the gap equation at fixed value of q. In the plot, Sfi is fixed to the numerical value O.TAq with Aq = 80 
MeV; but changing Sfi does not change the picture qualitatively. For q of the same order of Aq, the shape of the 
gap function is very close to that of the strip. In this case, the largest contribution to the gap comes from the lowest 
order harmonic n = 1 in Eq. ()112p . However, for small values of g/Ao, the solution of the gap equation has a solitonic 
shape. This is evident in the case q = O.IAq, in which A{z) « ±Ao for one half-period, then suddenly changes its 
sign in a narrow interval. In this case, the higher order harmonics play a relevant role in the gap function profile. 

In the right panel of Fig. |8l we report the plot of the physical value of q as a function of Sfi. In the window in 
which q — 0, the ground state is the homogeneous BCS state. At the critical value Sfi = Sfic ~ 0.695Ao, the ground 
state has q 0. This signals the transition to the inhomogeneous phase. Since q can be arbitrary small in proximity 
of the transition point, from the left panel of Fig. [8] we read that the ground state will consist of a solitonic structure 
for Sfi sa S^c- As Sfi is increased, the physical value of q increases as well. Therefore, again from the left panel of the 
figure, we read that the solitonic structure will continuously evolve into the strip structure. 

In Fig. m taken from iNickel and Buballal (|2009[ ). we report the plot of the BCS free-energy (solid line) and of the 
inhomogeneous phase (dashed line) as a function of Sfi. The dotted line corresponds to the free energy of the FF 
phase. The baseline corresponds to the free-energy of the normal phase. With increasing mismatch a second order 
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FIG. 9 Free-energy density as a function of for the BCS phase (solid li ne), FF p hase (dotted line) and solitonic phase 
(dashed line). Free energies are measured respect to the normal phase. From lNickel and Buballa (200^ '). 
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phase transition takes plaee from the BCS to the sohtonic structure at a value (5/^/Abcs — 0.7 which is smaller than 
the value at which the transition from the BCS phase to the FF phase takes place. On the other hand, the transition 
from the inhomogeneous phase to the normal one takes place almost simultaneously to the transition from the FF 
phase to the normal one. 

In the case of the FF phase, the transition to the BCS state is first order and to the normal phase is second order. 
On the other hand, the transition from the soliton lattice to the BCS phase is of the second order. This is possible 
because the gap function (jll2|) naturally interpolates between the homogeneous state, corresponding to a single soliton 
with infinite period, and the cosine-shaped solution. 

Moreover, the transition to the inhomogeneous state (11121) to the normal pha se is found to be of the first order, in 



agreement with the results of the smearing procedure (ICasalbuoni et aLl . |2004 ). At this transition point, the profile 



of the gap function is of the cosine type, i. e. the strip profile discussed above. 

The results of this analysis are in disagreement with those of the GL expansion ( Bowers and Rajagopaj 2002 : 



Larkin and Ovchinnikovlll964l ). which predicts that the transition from the strip to the normal phase is of the second 



order. Since the GL expansion is expected to be exact in proximity of a sec ond order phase transition , it is important to 
understan d the origin of the discrepa ncy among the GL result and that of lNickel and Buballa To this end, the 

authors of Nickel and Buballa ( 20091 ) compute the thermodynamic potential, Q, for a gap function A(z) = Acos(2gz) 
with a fixed value of q, as it is customary in the GL studies. Their analysis reveals that in this case, beside the local 
minimum of located at small values of A, and which is captured by the GL expansion, a global minimum appears 
for larger values of A. Then the authors argue that in order to capture this true minimum, which is responsible for 
the first-order transition to the normal phase, terms of at leas t eight h order in the GL expansion, which are usually 
neglected, should be included. In a subsequent paper (jNickei \2()m Nickel discusses the inhomogeneous phases for 
lower dimensional modulations in the NJL model and in the quark-meson model. His results confirm the replacement 
of the first order transition of the phase diagram of the homogeneous NJL phase by two transition lines of second order. 
These lines are the borders of an inhomogeneous p hase and they intersect a t the critical point. An interesting point is 
also the relation to the chiral Gross- Neveu model ([Gross and Neveul . 119741 ) . A more comp lete program would require 
inclusion in the study of higher dimensional modulations of the i nhom ogeneity (as done in lAbuki et al. ( 2012 ) for the 
chiral condensates), and, as pointed out by iBuballa and Nickell ( 2010l ). simultaneous study of both superconducting 
and chiral condensates. 



H. An aside: Condensed matter systems 



Fermionic systems consisting of two different "flavors" with mismatched Fermi surfaces and non-homogeneous 
condensates are quite generic and appear in various contexts. Particularly interesting are the population imbalanced 
superfluid (or superconducting) systems that can be realized and studied in laboratory. Examples of these type are 
gases of cold atoms, type-II cuprates and organic superconductors. The study of these systems allows to shed light 
on various aspects of superfluidity of asymmetric systems in a framework under experimental control. 

Of considerable importance for their similarity to quark matter are ultracold systems co nsisting of ferm i ons o f two 
different species, '0i and 'i/'2, c orresponding to two hyperfine states of a fermionic atom, see iGiorgini et ali (j2008( ) and 
Ketterle and ZwierleinI (|2008l ) for reviews. These fermions have opposite spin and one can change the number of up 
and down fermions at will. One of the most intrig uing aspects of the se systems is that the interaction between fermions 



can be tuned by employing a Feshbach resonance ( Chin et al . 2010l ). and therefore the crossover between the BCS and 



the BEG superfluid phases can be studied. This research field has grown in an impressive way in the last two decades 
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A mismatch between the populations of electrons in a superconductor can also be produced by Zeeman splitting. 
The magnetic field that couples with the spins of the electrons can be an external one or an exchange field. However, 
an external magnetic field couples with the orbital motion of the electrons as well, destroying superconductivity or 
leading to the creation of a vortex lattice structure. 

In order to reduce the orbital effect one employs 2D supercon ductors, i.e . filrn s of superconducting material or 
systems with a layered structure, and an in-plane magnetic field ( Bulaevskil. 1973). Go od candidates for LOFF su- 
perconductors of this type are heav y- fermion compounds like CeRu ^ |Huxlev et all . 1993 ): recently interesting results 
have been obtained with CcCOINs (jMatsuda and Shimaharall2007l) . Quasi- two-dimensional organic superconductors 
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(jUii et ali . 120061 ) like k-(ET) or A-(ET) salts, in particular A-(BETS)2FeCl4, are promising candidates for realizing 
the LOFF phase. High Tc superconductors are good candidates as well. In case one uses 3D materials one has to 
understand how the vortex lattice structure changes in presence of a LOFF periodic order parameter. 

As a final rem ark, we recall that the gapless CSC phases are the QCD analogue of a conden s ed matter pha se, known 
as Sarma ph ase (iForbes McNeil et a/.l. l2005l: ICubankova et aLl. 120031; iLiu and Wilczekl. 120031: ISarmal. Il963l). found to 



Wu and Yid . l2003h in the weak coupling limit 



be unstable (ICubankova et ali l2010l l200d iMannarelli et aLl . l2006at |Pao et all l2006l : ISheehv and RadzihovsWi l2006t 



III. THE THREE-FLAVOR NON-HOMOGENOUS PHASES 

In the previous chapter we have discussed the crystalline non-homogenous phases that can be realized in two-flavor 
quark matter. A plethora of possible structures have been analyzed by a Ginzburg-Landau expansion and by the 
so-called smearing procedure in order to determine the favored thermodynamic state. 

Here we extend the analysis to the three-flavor case which is relevant if the effective strange quark mass, Ms, is 
not too heavy. In compact stars Ms, will lie somewhere between its current mass of order 100 MeV and its vacuum 
constituent mass of order 500 MeV, and therefore is of the order of the quark number chemical potential /i, which 
is expected to be in the range (400 — 500) MeV. Furthermore, deconfined quark matter, if present in compact stars, 
must be in weak equilibrium and must be electrically and color ncutraH. All these factors work to separate the Fermi 
momenta of the quarks and thus disfavor the cross-species ECS pairing. 

In order to understand how the mismatch among Fermi momenta is linked to /3— equilibrium and neutrality, let us 
consider quark matter composed by u, d and s quarks with no strong interactions. We treat the strange quark mass 
as a parameter and we consider sufficiently long time-scales for which weak equilibrium and electrical neutrality are 
relevant. 

Since we are assuming that color interactions are absent, the chemical potential of quarks is diagonal in the color 
indices and the free-energy of the system is given by 

3 f^" 3 f^^ 

5^unpaircd(/i, ^J■e, Mg) = — / (P ~ tl-u)p^ dp + — (p - ^^d)p^dp 

Jo ""Jo 
^2 r {V^TMl-ps)p'dp+^ f'^\p-fi,)p'dp, (113) 







where the Fermi momenta are given by 

Pu=Pu P! = Pd Pf = Vf^'s - . (114) 
Weak equilibrium relates the chemical potentials of quarks with different flavors as follows 

2 11 , ^ 

Pu^p-^fJ-e ^ld = ^^.+ ^^^e /i5=/i+-/Xe, (115) 

and the electrical chemical potential is obtained from the neutrality constraint 

1^=0. (116) 

Solving Eq. (jll6l) one determines pe and by Eqs. (|114p one obtains the values of , and Pf for the electrically 
neutral unpaired phase. The effect of Ms is to reduce the number of strange quarks, which must be compensated by 
a larger number of down quarks to ensure electrical neutrality. It follows that there is a hierarchy of Fermi momenta 

Ps^<P^<P[, (117) 

and in Fig.[TU]we report the difference of the Fermi momenta Pf — P^ and P^ — P^ as a function of Ms, for /i ~ 500 
MeV. 

For Ms ^ p., the effect of a nonzero strange quark mass can be taken into account by treating the strange quark 
as massless, but with a chemical potential that is lowered by M^/(2/i) from /i -|- ^e/3. To the order Ms / p electric 



^ Actually, as already discussed in SecU] quark matter must be in a color singlet, however it has been shown in lAmore et al. 1 1 12002!) that 
projecting out color singlet states into color neutral states does not lead to a large change of the free energy. 
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FIG. 10 Fermi momenta dilferences P£ — , black solid line, and Pjf — Pf, red dashed line, as a function of Ms, for /i = 500 
MeV. 



neutrality, Eq. (|116|) . requires that fie — In this case P£ — ~ — Pf , and we need no longer to be careful 
about the distinction between Pp and /x, as we can simply think of the three flavors of quarks as if they have chemical 
potentials 



fJ-d = t-iu + 2 Sfi 



Pf 



/.is = /iu - 2 (5^ , 



(118) 



with 



M2 



(119) 



and we can write the chemical potential matrix as 



/^U,a/3 = Sal3 "S) diag (^„, fid, fis)i 



(120) 



The above results are strictly valid for unpaired matter, but one expects that, if the pairing interaction channel is 
not exceedingly strong, then it should still be qualitatively correct. In particular, it should happen that the hierarchy 
of Fermi momenta is the one reported in Eq. (|117p and the splittings those reported in FigfTUl 

As discussed in the two- flavor case, a small values of the chemical potential difference (5/z, i.e. of Ms, cannot 
disrupt the BCS pairing and the CFL phase will be energetically favored. However, when Mg is sufficiently large, 
the mismatch between the Fermi momenta becomes large disfavoring CFL pairing, that must be superseded by a less 
symmetric phase. 

As shown in Fig. [TUl for values of Mg comparable with fx the strange quarks decouple and only the 2SC phase wit h 
the pairing between up and down quarks seems realizable. However, as pointed out in lAlford and Raiagopall (120021 ) . 
once the constraints of electrical and color neutrality and /3-equilibrium arc imposed, the 2SC phase turns out to be 
strong ly disfavored or even excluded at least at zero-t emper ature (for finite temperature evaluations see for example 
Abukii (|2003h '). The authors of lAlford and Raiagopall (|2002[ ). performing an expansion in terms of the strange quark 



mass (at the leading non-trivial order), found that whenever the 2SC phase is more favored than unpaired matter, then 
the CFL is even more favored. The point is that the CFL phas e is extremely robust, because it allows pairing between 
all the three-flavor species. This analysis has been redone in Steiner et al\ (|2002l ) evaluating the density-dependent 
strange quark mass self-consistently. The results of lSteiner et al\ ( 2002 ) almost confirm previous conclusions, finding 
the 2SC favored in a very narrow region of density, that is likely to disappear once the hadronic phase boundary is 
properly taken into account. Notice that these results are valid in the weak-copuling approxima tion. Computation 
using NJL-like models with a stronger coupling (jAbuki and Kunihird . 120061 : iRuester et a/.l . l2006al ). see also Sec. IH.Bl 
find that a large portion of the phase diagram is occupied by the 2SC phase. 
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FIG. 11 (Color online). Left panel: gap parameters A3, A2, and Ai as a function of M^/fi for fi — 500 MeV, in a model 
where Ao = 25 MeV (see text). Right panel: free-energy of the CFL/gCFL phase, relative to that of neutral noninteracting 
quark matter and that of the 2SC/g2SC and 2SCus phases. At M^/fj, « 47.1 MeV (vertical dotted line) there is a continuous 
phase transition between the CFL phase and the gCFL phase. A bove /fi ~ 130 M eV (vertical dash-dotted line) unpaired 
quark matter has a lower free-energy than the gCFL phase. From lAlford et all (|2005d ). 



A. The gapless CFL phase 

Given that the effect of non-zero strange quark mass and electrical neutrality constraints is to pull apart the Fermi 
spheres of different quark species, there is little motivation to assume the symmetric CFL pairing reported in Eq. ([2|). 
Since the gap parameter gives the strength of the interaction channel, one might expect that a mismatch between 
the Fermi momenta reflects in a reduction of the interaction channel and thus in a reduction of the corresponding 
gap parameter. Therefore, the gap parameters should be now flavor dependent and one can consider the generalized 
pairing ansatz 



(121) 



1=1 



where Ai , A 2 and A3, describe down-strange, up-strange and up-down Cooper pairing, respectively. In lAlford et al\ 



(|2004 |2005g) the sup erconducting phase char acterized by the ansatz (|12ip has been studied by a NJL-like model at 
zero temperature (see iFukushima et al\ (j2005() for a study at nonvanishing temperature) , considering the in- medium 
strange quark mass as a free parameter, while the light quarks are taken massless. The gap equations, coupled to the 
neutrality conditions, have been solved and the corresponding free-energy has been determined. Besides the strange 
quark mass, the final results depend on the quark chemical potential, which is fixed to the numerical value n = 500 
MeV; moreover, the strength of the NJL-like interaction is fixed by the value of the homogenous CFL gap, Aq. The 
numerical value Aq = 25 MeV, corresponding to th e weak coupling regime, has been chosen. 

Some of the results of lAlford "oLI (|2004 120053) are summarized in Fig. 1111 On the left panel of the figure, the 
solutions of the gap equations as a function of / fi are plotted. As it can be inferred from the figure, at / ^ 47.1 
MeV, there is a continuous phase transition between the CFL phase and a phase characterized by A3 > A2 > Ai > 0. 
Note that the hierarchy of the gap parameters is in agreement with the hierarchy of the Fermi momenta splitting 
reported in Fig. I10( which substantiates the reasoning presented above about the relation between mismatched Fermi 
momenta and gap parameters. Actually, the mismatches of chemical potential are not only due to the finite value 
of the strange quark mass, but as well as to the nonzero values of the color chemical potentials fi^ and /ig, which 
however do not overturn the hierarchy in Fig. 1101 

The phase appearing at M"^ / fjL k, 47.1 MeV is dubbed gapless CFL, or gCFL, phase, and is the three- fiavor analogous 
of the g2SC phase discussed in the previous chapter. The symmetry breaking pattern is the same of the CFL phase, 
see Eq. but the spectrum of s ome fer mionic excitations is gapless. The possibility of a gapless CFL superconductor 
was first argued in lAlford et al\ (|2000a[ ). where a toy-model was used to infer the effect of a heavy Mg on the quasi- 
particle spectrum and it was found that if the condensates of light and strange quarks are below a certain critical 
value there is not a minimum excitation energy in the quasiparticle spectrum. As already discussed in the previous 
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chapter, the mechanism at the origin of gapless modes in a superconductor is quite general, and is due to a mismatch 
6fi between condensing fcrmions of the order of A. 

In case of three-flavor quark matter, the situation is comphcated by the presence of several interaction channels. 
In the gCFL phase the dispersion laws are more complicated than those reported in Eq. ([T0|) . because they depend 
not only on the strange quark mass, but on the nonzero values of the color chemical potentials fj,^ and fxs, as well. 
However, the result is rather similar, as gapless excitations appear in the spectrum. 

Comparing the free-energy of the gCFL with other candidate phases, e.g. the 2SC phase, or the 2SC+s, which is a 
three-flavor model in which it is taken into account the decoupled strange quark, or the so-called 2SCus in which only 
the u-s pairing survives {i.e. A2 > 0, Ai = A3 = 0) or the g2SC discussed in the previo us chapter, the gCFL phas e 



turns out to be energetically favored in a quite large window of the control parameter ( Alford et al. , I2OO 1 l2005ct) 



The comparison of the gCFL free-energy with the free-energy of so me of theses pha ses is reported in Fig. [TT] Besides 
the considered phases, other pairing patterns have been proposed (jlida et corresp onding to the uSC phase 



(with Al = and_A2, A3 ^ 0) and the dSC phase (with A2 = and Ai, A3 ^ 0), see (jFukushima et all 12005 



Ruester et all l2006al) for a comparison of the thermodynamic potentials of the various phases. 



However, as discussed in the previous chapter, the fact that a phase is energetically favored with respect to other 
phases is not a sufficient condition to ensure its stability. Indeed, the analysis of fluctuations around the mean fleld 
solutions may reveal that the phase does not correspond to a minimum of the free energy. As in the case of the g2SC 
phase, the gCFL phase is indeed chromo-magne tically unstable, because four of the e i ght M eissner masses of giuons 



become imaginary when gapless modes appear (jCasalbuoni et all l2005bl : iFukushimal . 120051) . Further increasing the 



mismatch among the Fermi surfaces, the masses of the remaining four giuons b ecome imag i nary, as well. The effects 
of temperature on the Meissner masses in the gCFL phase has been studied in IFukushimal (|2005|) . where it is found 
that for sufficiently high temperatures (of about 10 MeV for the parameter choice of Fukushima ( 20051 )). the gCFL 
phase becomes stable. However, such a temperature is much larger than the typical temperature of compact stars and 
therefore cannot be used as an argument in favor of the gCFL phase in these systems. Also in this case the instability 
can be associated to the tendency of the system to develop supercurrents, that is to realize a FF-like state. Then, one 
generalizes the FF state to more complicated structures trying to single out the correct ground state among various 
crystalline phases. We shall investigate various crystalline structures in the following sections. 



B. Three-flavor crystalline phase: Two plane waves 



In this section we present the evaluation of the simplest crystalline superconducting phase using a modified Nambu- 



Gorkov form alism. The app roach presented is based on the High Density Effective Theory (HDET) of lCasalbuoni et al. 



(|2001al ) and iNarduUil (|2002[ ) and the evaluation is performed with and without making a GL expansion. As already 
discussed in the previous chapter, the GL expansion of the free-energy in powers of the gap parameter A is under 
control close to a second order transition, so that we should consider only small value of the gap, or more precisely 
A/Aq ~ A/S^ — > 0. This assumption implies that pairing does not significantly change any number densities. In this 
context neutrality is achieved with the same chemical potential /ig = Af^/(4/i) and /i3 = = as in unpaired quark 
matter. We will comment briefly in a next subsection on the latter approximation. 



1. Nambu-Gorkov and HDET formalisms 

In the framework of the Nambu-Gorkov formalism one defines an multi-component Nambu-Gorkov spinor 



X(P)=I 7t . ) ' (122) 



such that in this basis the pairing between fcrmions with momentum p and —p is described by an off-diagonal term 
in the fermion propagator, see e.g. Eq. (|20|) . Therefore in the Nambu-Gorkov space, the quark propagator has a 
compact form. In the two-fiavor case, considering the FF ansatz (|4ip . we are considering a condensate of the form 

(V'«(^)C75V'd(2^)) (X Ae2*'^■^ (123) 

where we have for simplicity suppressed color indices. In this case we are considering pairing between u quarks with 
momentum p + q and d quarks with momentum p -\- q. This could be described via a propagator with terms in 
it that are off-diagonal in momentum space, rather than merely off-diagonal in Nambu-Gorkov space. However, it 
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greatl y simplifies the calculation changing to a basis in which the Nambu-Gorkov spinor is written as in lBowers et al. 
(l200lh 



x{p) 



(124) 



Mp-q) 

i'li-p-q) 

\^Ji-p + q) / 

The pair condensate is then described by terms in the fermion propagator that are ofF-diagonal in Nambu-Gorkov 
space, occurring in the ipu-'ipj and ipd-'tpu entries. In this basis the fermions that pair have opposite momentum, 
making the propagator diagonal in p-space and the calculation tractable. 

One must always keep in mind that the fermions that pair have momenta p + q and —p + q and not p and —p. Since 
p is an integration variable, we can shift it by p — >■ p — q thus rewriting the Nambu-Gorkov spinor in the two-flavor 
phase as 



x(p) = 



f tpuip) 

iy{p - 2*7) 

\^Ii-P + 2q)J 



(125) 



Once we turn to the three-flavor case, the simplest form of inhomogcneous pairing corresponds to the condensate 



A 



3 

E 

1=1 



(126) 



meaning that for each pairing channel we assume a Fulde-Ferrcll ansatz. with 2qj representing the momentum of the 
pair. 

In this case it is not possible to diagonalize the propagator in p-space, because three different fields are locked 
together in pairs and thus one can only eliminate in the off-diagonal terms of the propagator two independent momenta. 
However, as shown in (jllSp and pi9p . the separation between the s and d Fermi spheres is twice the separation between 
the d and s Fermi spheres and the u and s Fermi spheres, thus it is reasonable to expect Ai <C A2, A3. As a first 
approximation one can consider Ai = and in this case it is possible to diagonalize the propagator in momentum 
space. 

The form of the two-flavor Nambu-Gorkov spinor (|125p immediately suggests that we analyze the three-flavor 
crystalline phase with condensate (jl26p with Ai set to zero by introducing the Nambu-Gorkov spinor 



x(p) 



/ Mp) \ 

■0d(p- 293) 
■0s(p- 292) 
^^Ui-P) 

^J(-p + 2q3) 

vv;n-p+2q2)/ 



(127) 



It is the clear that it would not be possible to use this method of calculation if Ai were kept nonzero, except for 
the special case in which qi ~ q2 — qs- (That is, except in this special case which is far from sufficiently generic, it 
will not be possible to choose a Nambu-Gorkov basis such that one obtains a propagator that is diagonal in some 
momentum variable p.) Moreover, it seems unlikely that this method can be employed to analyze more complicated 
crystal structures. Indeed, it seems that the pairing with Ai = and A2 and A3 each multiplying a single plane 
wave, is the most complex example that is currently know how to analyze without usin g the GL expan sion. 

We now implement the calculation in the basis (jl27l) using the HDET formalism (INarduUi 120021) . We Fourier 
decompose the fermionic fields as follows: 



dfl 



(128) 



where i? is a unit three-vector whose direction is integrated over, fc, are three vectors, one for each flavor, that we 
shall specify below and ijj'j'^^x) (resp. ip^^ix)) are positive (resp. negative) energy projections of the fermionic fields 
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with flavor i and color a indices, as defined in ICasalbuoni et ali (|2001al ) and iNardullil (|2002l ). Note that these fields 
are written in a mixed notation, meaning that they depend both on the space coordinates and the unit vector v, 
which points to a particular direction in momentum space. Indeed, in the standard HDET approximation (employed 
to describe the homogeneous phases), fc^ = and the field ipfvi^) used to describe a quark in a patch in momentum 
space in the vicinity of momentum P = fiv. The introduction of the ki vectors means that now ipTvi^) describes a 
quarks with momentum in a patch in the vicinity of momentum fiv + ki and the chemical potential differences are 
then given by 



S^ii{v) = /ij - ^ - ki • V . 



(129) 



To reproduce (|127|) . then we choose 



ku^O 
kd = 2^3 

ks = 2q2 . (130) 
At the leading order in l//i (i.e. neglecting the contribution of antiparticles) the QCD Lagrangian can be written as 



47r 



(131) 



where 5^i{v) = Pf — fi — ki ■ v and the four vectors and (the latter used only below) are defined by = {l,v) 
and V" = (1,-v). 

To take into account diquark con densation, we add to (|131|) a NJL-like interaction term, which is then treated in 
the mean field approximation (as in Afford et al. ( 2005cl )). giving 



^^--IJ ^Kfi^f-.C^.i^le-'^^'^'+'^^y^ +h.c. , 



(132) 



where the condensate is given in (jl26p . It is convenient to work in a new basis for the spinor fields defined by 



where the unitary matrices Fa are given by 



^^1 

-^4,5 



1 

V3 



A=l 



F2 — — ^3 



1 



6,7 



Ta ± iTr, 



V3^ 



(133) 



V3 



with Ta = Xa/'2 the SU(3) generators and Iq the identity matrix. Introducing the velocity dependent Nambu-Gorkov 
fields 



V2 1 cr-. 



(134) 



and replacing any matrix M in color-flavor space with A^ab = Tr[i^jMFg], the NJL Lagrangian density can be 
written in the compact form 



A-u 1 



(135) 



where the inverse propagator is given by 



S 



AB 



{V ■ t +5^ia{v))5ab -Aab 

-Aab {V-£- dflA{-v))5AB 



(136) 



where = {£o,£\iv) is a four-vector. Here, i\i is the "radial" momentum component of ip, parallel to v. In HDET, 
the three-momentum of a fermion is written as (/j, -|- i\\)v, with the integration over momentum space separated into 
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an angular integration over v and a radial integration over —5 < £\\ < 5. Here, the cutoff S must be smaller than /i 
but must be much larger than Aq (and dfj,), see Fig[TJ 

If Ai = 0, the space dependence in the anomalous terms of the propagator can be eliminated by choosing 

fc-u + fed = 2q3 

k^ + ks = 2q2 , (137) 
and the calculation is technically simplified. It has been numerically checked in Mannarelli et al. ( 2006b[ ) that different 



choices of fc„, kd and fc^ satisfying (|137p yield the same results for the gap parameter and free energy. The different 
choices yield quite different intermediate stages to the calculation; the fact that the final results are the same is a 
nontrivial check of the numerics. 



From the Lagrangian (|135p . following a derivation analogous to that in lAlford et al\ (j2005a ). the thermodynamic 
potential per unit volume can be evaluated to be 



'47r2 



18 

E 

a=l 



+S 



din 



47r 



2A2 



Me 



127r2 



(138) 



where we have set A2 = A3 = A and where G is the dimensionful coupling constant. The dependence on G and of 
the ultraviolet cutoff 5 can be eliminated by using the CFL gap equation, through the parameter Aq defined by 



Ao = 2^/3 5 cxp 



2G^i^ 



(139) 



where Aq is the CFL gap parameter for Mg ~ and /Xg = 0. In Eq. (|138p . the Ea are the energies of the fermionic 
quasiparticles in this phase, which are given by the 18 roots of det5"^ = 0, given in Eq. (|136|) . The doubling of 
degrees of freedom in the Nambu-Gorkov formalism means that the 18 roots come in pairs whose energies are related 
by £'a('W,^||) = Eb{^v, £\\). One set of nine roots describes (V'li.ci, '0-i',«) ^-nd (^i;,s, pairing, while the other 

set describes {ipv,u,''J-'-v,d) and {ip^^u,4'-v,s) pairing. Since v is integrated over, the free-energy can be evaluated by 
doing the sum in (|138p over either set of nine roots, instead of over all 18, and multiplying the sum by two. 

The lowest free-energy state is determined minimizing the free-energy given in Eq. (|138p with respect to the gap 
parameter A and with respect to q2 and q^. One could also determine self-consistently the values of fie, and /ig 
which ensure electrical and color neutrality. Moreover, one could allow A2 =/= A3 and minimize with respect to the two 
gap parameters separately. However, in the results that we shall present in the next section we shall fix A2 = A3 = A, 
/ie = Mg/(4^) and /i3 = /ig = 0, as is correct for small A. 



2. Ginzburg-Landau analysis 



As discussed in the two-flavor case, the GL expansion allows to determine the gap parameter and the free-energy 
for a generic crystalline structure. We shall discuss the GL results for various crystalline phases in Section [111.01 in 
the present section we restrict to the case of the condensate reported in Eq. (|126p . We recall that the GL expansion 
is controlled by the ratio A/Aq ~ A/dfx and is reliable in the vicinity of a second order transition; this restriction 
implies that pairing does not significantly change any number density, and thus one can assume fj.3 = /ig = and 
lie m M"^ /AfjL as in the normal phase. 



For the condensate in p26|l . the GL expansion of the free-energy is given by (jCasalbuoni et aLl . l2005al : lMannarelli et al. 
2006bl) 



VL 




0(A6), 



(140) 



where f2„ is the free-energy of the normal phase and ar, Pr and Prj c a n be d etermined by the microscopic theory. 
Employing the HDET formalism, it has been shown in lCasalbuoni et oL ( 2005a) that ar and Pr hav e the same formal 
expression derived in the GL analysis of the two-flavor FF structure ( Bowers and Raiagopall 12002 ) , meaning that 



ai =aFF{qi,Sfii) 



,,2 



2qi 



loa 



qi + Sm 



qi - Sfii 



■log 



4(g|-<5M?) 



(2^Ao)2 



r2 qj 



2 ' 



(141) 
(142) 
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and in the weak coupling limit 



(5^*1 = -IJ,e , 



(143) 



The terms /3/,/, with I ^ J , are instead characteristic of the three-flavor case as they represent the interaction term 
between different condensates. We have that 



Pi 



dv 



1 



27r2 J At: {it - 2qi ■ v ~ 2(5jUi) {ie - 2q2 ■ v - 2(5^2) ' 



(144) 



and /3i3 is obtained from /3i2 by the exchange q2 — > <73 and /i, -H> and /323 from /3i2 by qi — >■ and /ig O 
These are the only terms in which a dependence of the free-energy from the relative orientation of the qj can arise at 
this order. 

One should fix the norms qi by the minimization procedure and finally look for the most energetically favored 
orientations of the three vectors qi. This is in general a complex task as it would require the simultaneous minimization 
of the free-energy as a function of the angles between the three different q/. What is usually done to circumvent this 
complication is to propose different definite structures for the q/, and choose among them the one with the lowest 
free energy. 

It has been found in ICasalbuoni et ali (l2005al ) that the energetically favored solution corresponds to Ai =0 and 
A2 ~ A3 with q2 parallel to q^. As to the norms qi, since we work in the GL approximation, it is possible to neglect 
the C'(A^) terms in the minimization of ft. Thus, one has that daj/dqj = 0, which is identical to the condition for two 



flavo rs and therefore gives the result qj = 7]\6iii\ where 77 = 1/zg ~ 1.2 (jAlford et all 120011 : iLarkin and Ovchinnikov 



1964 . where Zg was defined in Ea. (j85|) : considering Eq. (jl43| ) this implies that us and du quark pair with equal 
momenta 



The free-energy can now be written as 



Q = ml = ml 



1 



(145) 



(146) 



where 



A = A2 = A3 



(147) 



is the only independent gap parameter. Minimizing the free-energy expression with respect to A we obtain for values 
of Sji where a is negative, the solution 



and the free-energy at the minimum is given by 



/3 + /323 ' 



(148) 



- — - f^7i 



2(/3 + /323) ' 

which still depend on the angle, 0, see Fig. [TH between 52 and 53, by /323. 



(149) 



3. Testing the Ginzburg-Landau approximation 



We are now able to evaluate the gap parameter and the free-energy for the "crystalline" color superconducting 
phase 



A 



Q/3 



y2+e2''9"^-'-£"^3g,,3) 



(150) 



with and without the GL approximation ( Mannarelli et at . 2006b[) : here we have written q2 = qn2 and 53 = qn^. 
We report the results obtained for /i = 500 MeV and Ag = 25 MeV. The calculations are made varying Mg, but we 
plot quantities versus Af^ / jj, because the most important effect of nonzero Mg is the splitting between the d, u and s 
Fermi momenta given in Eq. (|118|) . and we report the results for four values of the angle between n2 and n^: (f) ^ 0, 
27r/3, Ttt/S and 3l7r/32. The lines correspond to the GL analysis, where we have plotted A and f2 — fi„ of Eqs. (|148p 



) 




M^/|^ [MeV] 




M t\ , I , I , I , I 

"flO 120 130 140 150 

M^/)X [MeV] 



FIG. 12 (Color online). Plot of A/Aq (left panel) and of the free energy relative to neutral non interacting quark matter (right 
panel) as a function of /fi for four values of the angle <j> between q2 and qs. The various hues correspond to the calculations 
done in the GL approximation whereas dots correspond to the NJL calculation, done without making a GL approximation. The 
full lines (green online) and circles correspond to (/> = 0, the dashed-dotted lines (magenta online) and diamonds correspond 
to (/) = 27r/3, the dashed lines (r ed online) and squares c orrespond to <f) = 7tt/8, the dotted lines (blue online) and triangles 
correspond to = 3l7r/32. From lMannareUi et a/I l|2006bl '). 



and (|149|) and using Eq. (|14ip to relate a to 6^ and hence to M^//x. The points correspond to the NJL calculation 
without GL, where J7 of Eq. (jl38p has been minimized with respect to A. 

The NJL calculation has a second order transition at M^//x ^ 151 MeV, sec Fig.[T3 corresponding to ~ 0.754Ao, 
for all values of the angle (p. This result is in agreement with the GL calculation, in which the location of the phase 
transition depends only on a, which is independent of (/). Near the phase transition, where A/Aq and hence A/Sfi 
are small, there is good agreement between the NJL calculation and the GL approximation, as expected. When the 
GL approximation breaks down, it does so conservatively, under- predicting both A and jfi — r2„|. Furthermore, even 
where the GL approximation has broken down quantitatively, it correctly predicts the qualitative feature that at all 
values of / jj, the most favorable crystal structure is that with </> = 0. Therefore, the GL approximation is useful as 
a qualitative guide even where it has broken down quantitatively. 

It is evident from Fig. [12] that the extent of the regime in which the GL approximation is quantitatively reliable 
and is strongly </)- dependent. In the best case, which it turns out is = 0, the results of the GL calculation are 
in good agreement with those of the full NJL calculation as long as A/Aq < 0.25, corresponding to A/6fi < 0.35. 
For larger 0, the GL approximation yields quantitatively reliable results only for much smaller A. For example, 
with (j) — 3l7r/32 the GL calculation gives re sults in quantitative agr eement with the full NJL calculation only for 
A/Ao < 0.04, corresponding to A/611 < 0.05 (|Mannarelh et al\ . \200m . 



The free-energy of the favored crystalline phase, corresponding to n2 = has a lower frcc-cncrgy of the CFL and 
gCFL phases in a window 128 McV < M"^ / ^ < 151 MeV. Since the GL approximation is reliable over the widest 
domain precisely for the case with n2 \\ M3, we note from Fig. [12] that the GL analysis provides good approximation 
to the results of the full NJL calculations. These results strongly motivates the study of more complicated crystalline 
structures employing the GL approximation. 

We can now get a deeper insight about the favored orientation of the qi from the analysis of /323- As explained 
above, this is the only term in the GL expansion that depends on the relative orientation between n2 and and can 
be rewritten as 



/323 = ^/(0), (151) 

where 

/((/,) = ^e I ^- (152) 

J An [le — ijv • 922 — l)(«e — yyt) • 713 + 1) 

can be evaluated numerically and the result is plotted in Fig ll3l showing that the minimum value occurs at (/> = 0, 
that is for rii = ^2, in agreement with the fact that the configuration that minimizes /((/>) does as well minimize the 
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FIG. 13 I{4>), defined in Eq. (|152p . versus tlie angle (f> between the wave vectors q2 and qs. 

free-energy, sec Eq. (|149p . We note that although !{(/)) is an increasing function of (j> it depends very slightly on the 
angle except close to the value (f) = tt, i.e. when the two vectors are antiparallel and diverges. 

Although the divergence of /323 at (/) = tt is not physically relevant, since it is not the value that minimizes the free- 
energy, it is worth considering this case to gain a qualitative understanding of the behavior of the GL approximation. 
We see in Fig. [14] that there are two pairing rings on the up quark Fermi surface, because some up quarks pair 
with down quarks forming Cooper pairs with wave vector 2^3 and other up quarks pair with strange quarks forming 
Cooper pairs with wave vector 2q2. However, as shown in the right panel of Fig. I14[ if = tt the two pairing rings on 
the up quark Fermi surface are close to be coincident. In the weak-coupling limit in which dfi/fi — ?■ (and Ao — > 
with (5/i/Ao fixed) these two rings become precisely coincident. We attribute the divergence in (523 to the fact that 
antiparallel wave vectors pay an infinite free-energy price and hence are forbidden, because of the coincidence of these 
two pairing rings. In contrast, if ~ 0, the two pairing rings on the up Fermi surface are as far apart as they can 
be, and P23 and the free-energy of the state are minimized. This qualitative understanding also highlights that it is 
only in the strict GL and weak coupling limits that the cost of choosing antiparallel wave vectors diverges. If A/Sfi 
is small but nonzero, the pairing regions are ribbons on the Fermi surfaces instead of lines. And, if is small but 



FIG. 14 (Color online). Sketch showing where on the Fermi surfaces pairing occurs for condensates in which q2 and qs are 
at an angle (f). The light gray (red online) ribbons on the d and u Fermi surfaces indicate those quarks that contribute the 
most to the (ud) condensate with gap parameter A3 and wave vector qa, which points upward in both panels. The dark gray 
(blue online) ribbons on the u and s Fermi surfaces indicate those quarks that contribute most to the (us) condensate with gap 
parameter A2 and wave vector q2. For descriptive reasons, in the left panel it is shown the case <j) = 70°, and in the right panel 
the antipodal case, with the two vectors pointing in opposite directions. We have greatly exaggerated the splitting between the 
Fermi surfaces, relative to the values used in the calculations reported in Sec. IIII.BI 
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not taken to zero (as of eourse is the case in Fig. [T4)) then the two ribbons on the up Fermi surface will have slightly 
different diameter, as the figure indicates. This means that we expect that if we do a calculation at small but nonzero 
Aq ~ Sfi, and do not make a GL expansion, we should find some free-energy penalty for choosing = tt, but not a 
divergent one. This is indeed the result that one obtains without using the GL expansion, see FiglT^l These results 
also explain why the breakdown of the GL expansion in FigH^] happens for very small values of A for </> « tt, while 
only for larger values of A for = 0. Indeed, increasing (j) from to tt, the /323 term in the GL expansion increases, 
and therefore the radius of convergence of the GL expansion decreases. 

A more quantitative study of the radius of convergence of the GL approximation would require evaluating (at least) 
the A^ terms, whose coefficients we shall generically call 7. Because we are working in the vicinity of a point where 
a = 0, the first estimator of the radius of convergence that we can construct comes by requiring 7A^ ^ (/? + /323)A^. 
Thus, the results of the comparison in Fig. [T2]are not conclusive on this point, but they certainly indicate that the 
radius of convergence in A decreases with increasing 0, and tends towards zero for — > tt. 



4. Chromo-magnetic stability of the three-flavor crystalline phase 



Given that the three-flavor crystalline phase in Eq. (|150p is thermodynamically favored with respect to the CFL and 
normal phases, it re mains to be proven th at it is chromo-magnetically stable. This issue has been discussed within 
the GL expansion in Ciminale et all ( 2006| ). where the Meissner masses of gluons have been evaluated and it has been 
shown that this phase is indeed chromo-magnetically stable. To take into account gluons, the HDET Lagrangian has 
to be extended to the next-leading order (NLO) in obtaining in momentum space: 



f f 


a/3\ 


y ■e + 2ij_ 





where 



and 



(153) 



(154) 



(155) 



In the HDET there are two effective vertices describing the coupling between gluons and quarks. Either one gluon 
couples to a quark and an antiquark (three-body vertex, coupling ~ g) or two quarks couple to two gluons (four-body 
vertex, coupling ^ g^). These two vertices originate from the terms in Eq. (jl3ip with one and two momenta £. 




FIG. 15 (Color online). Squared Meissner masses of the gluons As (solid blue line), (dotted green line). As (dashed black 
line), and Ai (dot-dashed red line), in units of the CFL sq uared Meissner mass, vs M^/fi. Left panel: Longitudinal masses. 
Right panel: Transverse masses. Adapted from lCiminale et al., (,200^ 1. 
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Theses two vertices lead, at the order of , to two distinct contribution to the Meissner mass of gluons. The 
four-body coupling gives rise to the contribution / {2tt'^) , which is the same for all the eight gluons and is the 
same one has in the CFL phase. 

The three-body coupling gives rise to the polarization tensor: 



iWZix, y)= - Tr[z y) i H^i i S{y, x) z H'^] 



where the trace is over all the internal indexes; S{x,y) is the quark propagator, and 

\ 



(156) 



(157) 



is the vertex matrix in the HDET formalism. 

The Meissner masses are defined in terms of the gluon self-energy in momentum space, see Eq. ([57)) . and in this 
case the mass matrix of gluons in the adjoint sector a = 3, 8 is not diagonal. Therefore we introduce the two linear 
combinations 




Ai3 = cos O^Ais + sinOiAis , Ais = - sinO^Ais + cos O^Ais 



(158) 



which are eigenstates of the mass matrix. In (jl58p the subscript i denotes the spatial component of the gluon field; 
it is easily shown that the mixing angle satisfies the equation 



tan 26*, 



2Ml 



Ml., 



Ml 



the corresponding Meissner masses are the eigenvalues of the matrix 



(159) 



(160) 



which turn out to be positive and are reported in Fig. 1151 

On the left panel of Fig. [T2] we report the longitudinal (i. e . zz) components of the squ ared Meissner masses against 
M'^/fj,, in units of the CFL squared mass (iRischkel . l2000al: ISon and Stephanovl . l2000bl ). at the ©(A^); on the right 

panel the results for the transverse (i.e. xx) squared Meissner masses are given. 

These results arc the analogous of those obtained in the two-flavor case in Giannakis and Re3 (HooHE). Moreover, 
the transverse mass of Ag, although positive, is almost zero, being three order of magnitude smaller than the other 
ones. The conclusion is that the three-flavor condensate in E^ 
within the GL approximation 



The results of Ciminale et al. 



1501) has no chromo-magnctic instability, at least 
2006^ have not been extended to more complicated 



structures, however, for more complicated crystalline structures ( Raiagopal and Sharma . 2006[ ) the Meissner tensor 
should be positive definite for small values of A, since it is additive with respect to different terms of order A^ in 
the GL expansion ([Giannakis and Ren . l2005bl) . These considerations suggest that crystalline structures can be good 
candidates in removing the chromo-magnetic instability of the homogeneous gaplcss superconducting phases of QCD. 



5. Influence of 0(1/^) corrections 



Since the Meissner masses have been determined including the 0{1/^) corrections in the HDET Lagrangian, for 
a consistent calculation it is necess ary considering the effect of 0{l/fj,) corrections on the gap parameters and the 



free-energy (jCasalbuoni et all l200q ) . These corrections amount to a shift of the strange quark Fermi momentum to 



lower values with respect to the corresponding chemical potential 



Ps 



2fis 2fi \ 2/i 



(161) 



and thus increasing the difference between the u and s strange chemical potential, without affecting the d chemical 
potential. Thus (neglecting corrections proportional to /is and /ig which are of order 0{A^)) one has that 
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FIG. 16 (Color online). Left panel: Quark chemical potential difference as a function of //i. The dashed blue line corresponds 
to 5^2 and the dash-dotted red line corresponds to SjJis; their common value M^/(8/i), obtained neglecting 0{1/^) corrections, 
corresponds to the solid black line. Right panel: gap parameters as a function of M^/fi. The solid black line represents the 
solution A — A2 = A3 obtained neglecting 0(1/ f^) corrections; the dashed blue line and the dash-dotted red line represent 
respectively A2 and A3 obtained including 0(1/ ^i) correcti ons. All the gap parameters are normalized to the CFL homogeneous 
gap Ao = 25 MeV. Adapted from iCasalbuoni et al\ (|2006l '). 



= f , (162) 



and 6112 > (neglecting higher order effects, as in ICasalbuoni et al\ ()2005a|) . one would get the result 6^2 = Sfi^ 
since in this limit fi^ = M'^/A^). The results for the splitting of chemical potentials are reported in the left panel of 
Fig. [161 together with the result obtained in the large n limit. 

As a consequence of the 0(1/ fj.) corrections, the two gap parameters are not equal and A2 < A3, as shown in the 
right panel of Fig. [161 This effect is akin to the one considered in the gCFL phase, see Fig. [TTl but in that case the 
chromo-magnetic instability prevented the splitting of the gap parameters. The effect of the corrections considered 
here is to enlarge the LOFF window, as it can be seen from the right panel of Fig. [1^1 Neglecting 0(1/ fi) corrections 
the crystalline phase was characterized by one gap, and by condensation in two channels: u — s and u — d; including 
0(1/11) corrections there is pairing in both channels only for small values of M^/^, but the crystalline phase extends 
to larger values of / fi where only A3 7^ and therefore there is no u — s pairing. This phase is the analogous of 
the two-flavor FF phase, with a chemical potential difference determined by the ratio of the strange quark mass to 
the average quark chemical potential. 



C. Ginzburg-Landau analysis of crystalline structures 



As in the two-flavor case, more c omplicated crystalline structure can be considered and analyzed by the GL expan- 
sion (|RaiagODal and Sharmal 120061 ) . The most general crystalline diquark condensate was given in Eq. ([7]) , which we 
report here for convenience 



(oiv-rL^fjo) « 



3 

E 

1=1 



(163) 



This condensate is antisymmetric in color (a,/3), spin, flavor (i,j) indices, and corresponds to a inhomogeneous 
generalization of the CFL ansatz, but now each gap parameter has a periodic modulation in space corresponding to a 
crystalline structure, e.g. A2 is associated to a crystal due to u-s pairing which is described by the vectors q™, where 
m is the index which identifies the elements of the set {92}; meaning that the q™'s are the reciprocal vectors which 
define the crystal structure of the condensate. 

Self-consistent computations based on Eq. (|163|l are very complicat ed even within the GL calcula tion scheme; as a 
consequence, several assumptions and simplifications have been used in lRaiagopal and Sharma (Soo^. As discussed in 
the previous section, the pairing between d and s quarks can been neglected; this is a very reasonable approximation. 
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since the difference of chemical potentials between d and s quarks is approximately twice the imbalance of chemical 
potentials between u — s and u — d. Therefore, Ai = 0, and one can take A2 = A3, neglecting 0{l/ii) corrections. 
Another simplification is to consider only crystalline structures with wave vectors {q^} and {<7™} with equal modulus. 
This means that any phase is determined by one gap parameter, b y the modulus of the to t al pai r momentum, q, and 
by the unit vectors which determine the two crystalline structures (jRaiagopal and Sharmal . [2006h . The corresponding 
pairing ansatz simplifies to 



(OICl^^lIO) (X AJ] • (164) 

A further simplification of iRaiagopal and Sharnia ( 2006 ) is to consider only structures which arc exchange sym- 
metric, which means that {n2} and {^3} can be exchanged by some combination of rotations and rcficctions applied 
simultaneously to all the wave vectors. 

The thermodynamic potential within the GL expansion up to the sextic order in the gap parameters is given by 



f7(A2,A3)==^ 
2^2 



Pa(|A2| 



|A.3|2) 



r(|A2 



A,\^) + l{\A,f + \A,f) 



J23 



A2nA3p 



7233 



A2PIA3I4 



7322 



Ao|"|A3p 



(165) 



where we have explicitly diversified A2 and A3 in order to show the interaction between the condensates in the (ud) 
and the (us) channels . The coefficients a, (5 and 7 are the same coeflacients computed for the two-flavor case in 
Bowers and Raiagopall ( 2002 ). see the expression in Eq. (pi)) . the remaining terms are peculiar of the three- flavor 
case. In the numerical computations ri(A, A) is minimized to compute the physical value of A and for any crystalline 
structure considered it is found that the thermodynamic potential is always bounded from below. Therefore, it is 
possible to compute the numerical value of A and the value of the free-energy densi ty at the minimum. 

Among the many crystalline structures considered in IRaiagopal and Sharmal ( 20061 ). two of them have been found to 
have the lowest free-energy: they are called the CubeX and the 2Cube45z structures. In the CubeX crystal, {^2} and 
{n^} each contains four wave vectors, pointing to the vertices of a square; the eight vectors together point toward the 
vertices of a cube. In the 2Cube45z crystal, {n2} and {ris} each contains eight wave vectors, pointing to the corners of 
a cube; the two cubes are rotated by 45 degrees about the z— axis. The computation of the fr ee-energy density at the 
global minimum shows that the CubeX crystal is favorable over all the structures considered in lRaiagopal and Sharma 
(|2006h in the range 



K'P 

2.9Ao < — < 6.4Ao 



similarly for the 2Cube45z crystal it is found that 



6.4Ao < — < IO.4A0 



(166) 



(167) 



Putting together the above results, the crystalline superconducting phase is favored with respect to the homogeneous 
CFL, gCFL and unpaired phases in the range 



2.9A0 < — < IO.4A0 



(168) 



The above condition can be translated to a condition on ^ only if the constituent strange quark mass as a function 
of ^ is known. This topic will be discussed in the next section. 

We remark that from the quantitative point of view, the result in Eq. (|168p should be taken with care. The 
GL expa nsion is reliable only if the gap parameter A is small, compared to J/i. On the other hand, the numerical 
results of IRaiagopal and Sharmal ( 20061) show that within this expansion, the ratio A/Sii turns out to be of order 
one. This clearly signals that the results lie beyond the validity of the expansion itself. Moreover, the coefficients 
in the expansion in Eq. (|165p depend on the microsc opic model. Therefore, Eg. (11681) is certainly model dependent. 
Nevertheless, the qualitative picture which arises in IRaiagopal and Sharma |2006[ ) seems to be quite robust: the 



crystalline superconductor has lower free-energy with respect to the single plane wave state. This conclusion seems 
very reasonable and model independent: in fact, the phase space for pairing in the case of CubeX and 2Cube45z is 
larger than the one corresponding to the plane wave structure; this suggests that the free-energy gain for multiple- 
planc-wavcs pairings are larger than the ones obtained for the single plane wave state. 
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1. LOFF window in the QCD phase diagram 



It is important to identify, within a self-consistent computation, the range of temperature/chemical potential in 
which the crystalline phase is expected to be thermodynamically favored. In the previous sections we have been 
able to identify the LOFF window in terms of the dimensionless ratio M'^/fj?; in order to translate this window 
to a wind ow in the comput ation of the in-medium strange quark mass is necessary. This investigation has been 
started in llppolito et alii ( 2007 ) considering vanishing temperature. In that article, the in-medium quark masses are 
computed self-consistently within a NJL model. This model offers a quite simple theoretical framework to discuss 
simultaneously dynamical breaking of chiral symmetry, as well as color superconducting p airing. 

I n the simple case of two plane waves, pSOp . it was found in ICasalbuoni et~al. ( 2005a ): Mannarelli et al. ( 2006b ): 
and lRajagopal and Sharinal ( 2006[ ) that this state is energetically favored for 



where 



4.8OA0 < — < 7.56Ao 



Ao = 2^/3(Ao-,)exp(-^ 



(169) 



(170) 



i s the CFL gap in the chiral limit evaluated by the NJL model. At the reference point considered in llppolito et al. 



(|2007n . Ao 643 MeV, G(Ao) ~ 1.81/Ag have been determined and taking n = 50 MeV one has that Aq « 50 



MeV. Using the results for the constituent quark masses computed in llppolito et ali ( 20071 ). equation ()169p can be 
translated to a window for the existence of the LOFF phase in the quark number chemical potential, namely 



467 MeV < < 488 MeV 



(171) 



Therefore there exists a small but finite window in fi where the structure (jl50p . is favored in comparison with the 
normal state. 

As discussed in the previous section, the analysis of Raiagopal and Sharma (" 2C)06|) sh ows that more complicated 
crystalline structures can lead to a smaller free energy. In lRaiagopal and Sharmal (2006) it w as found that a CCSC 
phase exists in the interval reported in (jl68p which using the self-consistent treatment by llppolito et oil (|2007l ) 
translates in 



442 MeV < ^ < 515 MeV 



(172) 



This result is certainly model dependent, however it shows that the actual extension of the LOFF window might not 
be very large and therefore it may occupy only a fraction of the volume of hybrid compact stars with a quark core. 
We shall discuss more on this topic in the next chapter. 



D. Shear modulus and Nambu-Goldstone modes 



The three-flavor CCSC phase has several low energy excitations; besides gapless quasifermions, there are scalar 
modes which are related to t he spontaneous break ing of translational symmetry, the phonons, of the chiral symmetry. 



the pseudoscalar NB bosons (jAnglani et ali 120071 ) and to the breaking of U{1)b symmetry, the phonon. 



In view of the astrophysical applications of the crystalline phases that we shall discuss in the next chapter, the 
low energy effective actions of the low energy scalar modes of the LOFF pha se have to be computed. The ph onon 
Lagrangian for several two- flavor crystalline structures has been determi ned in Casalbuoni a^.l (|2002bl . l2001al ). and 
discussed in lII.D.3l while for the thre e-flavor case it has be en studied in Mannarelli et aZ.I ~ 2007 ). Regarding the NG 
bosons, we shall discuss the results of lAnglani et al. (2007) obtained for three-flavor crystalline phases. 



1. Phonons Effective Action and Shear IVIodulus 



According to the basic theory of elastic media (jLandau and Lifshit the elastic moduli are related to the 

potential energy cost of small deformations of the crystal. Therefor e, the evaluation of the shear modulus requires the 



knowledge of the low energy Lagrangian for the displacement fields (jMannarelli et a/.l . 120071 ). These displacements can 



be described by phonons that originate from the spontaneous breaking of rotational invariance by means of spatial 
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modulation of the gap parameter A. More in detail, phonons are small position and time dependent displacements 
of the condensate and, since the three condensates that characterize the crystalline phase can oscillate independently, 
we expect there to be three sets of displacement fields Ui{x). In the presence of the phonons, then. 



and the Lagrangian that includes fluctuation on the top of the mean field solution can be written as 



1_ / i^ + i/l A"(a;) 
2-^1 A"(a;) {i^-i/i)' 



1 



X 



+ ^tr((A^)-A«) 



(173) 



(174) 



where tr represents the trace over color, flavor and Dirac indices. 

To find the low energy effective action which describes the phonons and the gapless fermionic ex citations, one should 
integrate out the high energy fermion fields. The procedure is detailed in lMannarelli et al\ (|2007l ). where the effective 
action is derived starting from a NJL-like microscopic model. The final form of the effective action, of the phonon 
fields, is given by 



iS[u] ^ \og{Z[u]) 



(fx 



^ tr((A")^A")l+iTr„glog(5-i) 



16G 



(175) 



where Trng stands for the trace over the Nambu-Gorkov, color, flavor and Dirac indices and a further trace over a set 
of functions on space-time containing all energy modes. The full inverse propagator, is given by 



i^ + l^ A"(.t) 
A"(x) [li^-i/if 



(176) 



and it includes interactions of the phonon fields. 

For the crystal structures CubeX and 2Cube45z, the full inverse propagator cannot be inverted, so a GL expansion 
has been performed in order to obtain the effective action for the phonon field, first separating the full inverse 
propagator into the free part, Sq^ , and a part containing the condensate, E, as follows: S^^ = S^^ + E, where 











and 



S = 



A"(r) 
A"(r) 



(177) 



(178) 



Then, one can expand the term log(S' ^) that appears on the right-hand side of Eq. (|175p as 



Tr„g(log(5o-i + S)) = Tr„g(log5„-^) + Tr„g(5oS) - iTr„g(5oS)' 



(179) 



The Trng(logS'Q ^) term is related to the free-energy of unpaired (normal) quark matter. Furthermore, it is easily 
found that only even powers of (5*0^) contribute to the trace over Nambu-Gorkov indices. Expanding the effective 
action in powers of 0, and including the first non-trivial quadratic term, which is calculated to order A^, one obtains 
the effective action 



EE 



(M 



(180) 



where k = k2 — ki is the four momentum of the phonon and 



1 



ij^ + iT + h + if^3)iJ>-iT + h-i(^k) 
1 



(181) 
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where the trace is over Dirac indices. Integrating the expression above one obtains 

M - ^ / d^x Y.^'^Y. [^o("" ■ «/)^o« • ui) - {nT ■ d){nT ■ ui){nT ■ d){nT ■ uj)] , (182) 

where 

where Zq was defined in Eq. ([55|) . 

This is the low energy effective action for phonons in any crystalhne color superconducting phase, valid to second 
order in derivatives, to second order in the gap parameters A/'s and to second order in the phonon fields u/. Note that 
there are no terms in (|182p that "mix" the different uj{k), meaning that at this order the displacement of the various 
crystals can be considered separately. This follows from the fact that the Lagrangian conserves particle number for 
every flavor of quarks, which corresponds to symmetry under independent global phase rotations of quark fields of the 
three-flavors, meaning independent phase rotations of the three A/. The effective action should be invariant under 
these rotations and hence A/ can only occur in the combination A^A/. The mixing between different fluctuation can 
only appear in higher order terms in A/, e.g. by by terms like iJ.^\AiAj\^dujduj/SiJ?. 

The coefficients ki dete rmine the potential ener gy cost of a fluctuation and therefore they are related to the shear 



modulus. In particular in iMannarelli et al\ (|2007t ) it is found that for the the two favored structures, 2Cube45z and 



CubeX, the shear modulus is a 3 x 3 non diagonal matrix in coordinate space with entries proportional to 



16 , , 

t'CQM = YK. (184) 



Evaluating k yields 



2 



.CQM = 2.47MeV/fm3(^^^j (^^) (185) 
and considering 

350 MeV < ^ < 500 MeV . (186) 

and 

5 MeV < A < 25 MeV , (187) 

one has that 

0.47 MeV/fm^ < i^cqm < 24 MeV/fm^ . (188) 
The standard neutron star crust, which is a conventional crystal of positively charged ions imm ersed in a fluid of 



electrons (and, at sufficient de pth, a fiuid of neutrons) has a shear modulus approximately given by (jMannarelli et at. 



20071 : IStrohmaver et aLl . ll99lf ) 



0.092 keV/fm^ < i^nm < 23 keV/fm^ , (189) 

thus, the crystalline quark matter is more rigid than the conventional neutron star crust by at least a factor of 20-1000. 
Note that the three-flavor CCSC phase is also a superfluid, by picking a phase its order parameter does indeed break 
the quark- number U{\)b symmetry spontaneously. These results demonstrate that this superfluid phase of matter is 
at the same time a rigid solid and a superfluid. 



2. Goldstone modes 



The superfluid property of the CCSC phase derives from the existence of a massless NGB associated to the breaking 
of U(1)b- A c tually, the crystalline condensate breaks the same global (and local) symmetries of the CFL phase 
(jAlford et all 119990 ). leaving unbroken a global symmetry group: SU {?))c+l+r x ^2, see Therefore, there 
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are nine NGBs due to the spontaneous breaking of U{1)b and of the chiral symmetry, but only the U{1)b boson 
(iJ— phonon) is massless. Indeed, the pseudo-scalars associated to chiral symmetry have mass, because this group is 
explicitly broken by quark mass terms and by chemical potential differences. 

Alike the phono ns descr i bed in the previous section, the NGBs discussed here describe the fluctuations of the 
condensate, see e.g. lEguchil(ll976l) . The pseudo-NGBs are related to fluctuations in flavor space, while the iJ— phonon, 
(j), describes fluctuations in baryonic number. Here we shall focus on the ff— phonon and introduce it by means of 
the transformation ip — >■ Wip with U = exp|zq!)/ /|, where / is its de cay constant. Taking into account the unitary 
rotations, the HDET Lagrangian takes the form (jAnglani et all 120071 ) 



C 



dn ^ / {V■£ + 6^lA{v)) Sab -^*ba 

-^AB {V■l-5^JiA[~v))5AB 



An 



XB +{L^R) , 



(190) 



where all the quantities have been defined in Sec. IIII.B.ll see Ea. (|136p . except S^s, which is given by 

Eab = AJ(r)Tr[e/(FAC/^)^e/FsC/t] . 



(191) 



Expanding the Lagrangian in the the scalar field one derives three body and four body interaction vertices. Then, 
at the leading order in A, there are two contributions to the phonon self-energy, due to a tadpole-like and a 
charm-like Feynman diagrams. Technically the calculation is very similar to the one sketched in the previous section 
and the result is that the charm-like diagram contributes to the self-energy by 



1=2 qT 



cj){~m{k) pr(fco,fc) 



(192) 



with k = (fco, k) and 'P^{ka,k) corresponds to the HDET version of Eq. (|18ip . 



1 



5^ 



{V ■ i + Sn- qj' ■ v)[V ■ {£ + k)+Sn-q''P -v 
1 

{V ■l-5ii-q'p ■ v)[V ■ {l + k) - 5ii- qf 



the tadpole-like contribution is given by 



Si 



tad 



*-|E^? E/(0i vTiko^o,k 



0) 



(193) 



(194) 



We do not provide further details of the calculation which can be found in lAnglani et al\ (l2007|), and we only quote 
the final result valid at small momenta: 



^k) = ^(bi~k) [kilo - hk,V,,] <t){k) 



where 



1 



1=2 qT 



47r (5/i-gr-^ + «0+)' 

d£l ViVj 

{6fi-qf -v + iO+y 



(195) 



(196) 



(197) 



Specializing these results to the case A2 = A3 = A, jq^l = 1 93 1 = 9 (|Raiagopal and Sharmal . [200i), and requiring 
canonical normalization of the Lagrangian in Eq. (|195p . leads to 
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where P is the number of plane waves comprising each crystal. Upon substituting this result in Eq. (|197p . and 
specializing it to P = 1, we find Vij = [diag(0,0, where we have chosen q2 and along the positive z-axis. For 
the two cubic structures corresponding to the values P = 4 (CubeX) and P = 8 (2cube45z) we find Vij ~ Sij/i 



I.e. 



the velocity is isotropic and has the value 1 / -s/S. 

The extens i on of the above procedure to the octet of pseudo-NGBs is straightforward, and is presented in 
Anglani et al. ( 20071 ). Here we just recall that one of the important result is that the squared masses of the 



pseudo-NGBs are always positive, thus kaon condensation for the considered crystalline phases is excluded, at least 
at the order A^. 



IV. ASTROPHYSICS 



The artificial creation of the low temperature- high density conditions ap propriate for testin g the properties of color 



superconductors is one of the challenging aims of high-energy experiments (jKlahn et a/.l . 120121 ). To date, however, the 



unique laboratories which may produce these extreme conditions are CSOs. Note that even if the relevant densities 
and low temperatures conditions were reached in a terrestrial laboratory, the conditions realized in a CSO are different 
from those produced in an accelerator, because in the former case quark matter is long-lived, charge neutral and in 
/3-equilibrium. 

In the inner co re of a CSO, the extremely high density and low temperatures may favor the transition from nuclear 
to quark matter ( Bavm and Chin . 1976 : Collins and Perrv . 1975 : Ivanenko and Kurdgelaidz^ 1969i 1965 ) and in turn, 
to the color superconducting phase. Indeed, if matter is compressed at densities about a factor 5 larger than the 
density of an ordinary nucleus, a simple geome trical reasoning suggests that baryons are likely to lose their identity and 
dissolve into deconfined quarks (|Webej . ll999l) . In this case compact (hybrid) stars featuring quark cores surrounded 
by a nuclear mantle would exist and they could provide a unique window on the properties of QCD at high baryon 
densities. 

Whether matter in the interior of compact stars is compressed by gravity to such extrem e densities is still an 
open quest ion. Recent ast r onom ical observation of very massive CSO (jPemorest et al . 2010l ) seem to disfavor this 



possibility (jLogoteta et a^.l . 120121 ). but the results depend on the poorly known equation o f state (EoS) of matt er at 



high density , and the possibility that hybrid stars of about 2Mq have a quark matter core (|Alford et all l2005al) 



CCSC core (jlppolito et all 120081 ) cannot be excluded 



or a 



Assuming that deconfined matter is present, it should be in a color superconducting phase, because its critical 
temperature (in weak coupling) is given by Tc ~ 0.57A, and, as we have seen in the previous chapters we expect 
that A ^ 5 100 MeV. Lower values of the gap parameter, of the order of the keV up to few MeV, are realized 
in the spin-1 sinlge f lavor pairing phase (jAlford et all l2003l 1 19981 : iBuballa et 'ail . 120031 : ISchafeil l200d ISchmittl . l2005t 



Schmitt et all 120021 ) . In any case, for the greatest part of the compact star lifetime, the temperature is much lower 



than this critical temperature and the color superconducting state is thermodynamically favored. 

The presence of color superconducting phases in the interior of the CSOs can be probed by a number of astrophysical 
observables linked to the microscopic internal properties. Unfortunately, the available observational data do not allow 
us to infer in a unique way the internal structures of CSOs, but the investigation of the astrophysical signatures 
can help to pave the path connecting theoretical models to real observations. Searching for "the smoking gun" of 
deconfined matter is a hard task that should take into account several results about mass-radius relations, equations 
of state, thermodynamic observables, r-modes, glitches and magnetic fields. In the following we will present a brief 
state-of-the-art of some astrophysical implications of color superconductivity and in particular, of a crystalline color 
superconducting core. 



A. Gravitational Waves 

As reported in Section IlII.DI the crystalline color superconducting phase is extremely rigid, with a shear modulus 
larger than the one of standard neutron crust by at least a factor 20. The existence of such a rigid core within neutron 
stars may have a variety of observable consequences. If some agency (for example magnetic fields not aligned with the 
rotation axis) could maintain the rigid core in a shape that has a nonzero quadrupole moment, gravity waves would 
be emitted by a spinning CSO. The size of the distortion can be measured by the equatorial ellipticity 



(199) 
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FIG. 17 (Color online). Exclusion plots in the A-/i plane for the Crab (left panel) and the pulsar J2124-3358 (right panel) 
obtained assuming that these stars have mass M = 1.4 M©, radius 7? = 10 km, uniform mass distribution and are maximally 
strained. The area above each line is exclu ded by the correspondin g model. The rectangular box is the theoretically allowed 
region of u and A for the CCSC phase, see iMannarelli et al\ (|2007l ) . Left panel: The dashed lines correspond to the results of 
iLinI I2OO3) which considers the S3/S4 LIGO data of the Crab pulsar, for the case (Jmax = 10"'' (upper heavy red dashed line) 
and (Tmax = lO"'^ (lower light blue dashed line) and assuming that the whole star c onsis t s of m aximally strained crystalline 
color superconductin g matter. The solid lines have been determined using Eq. (7) of|Li3 (|2007l ). but with the S5 LIGO data 
(|Abbott et all I2OIOI ) and considering the case of (Tmax = 10 (upper heavy red solid line) and amax = 10 (lower light blue 
solid line). Right panel: Dependence of the ex clusion region on the size of the crystalline core for the pulsar J2124-3358 obtained 
usmg the S5 LIGO data (|Abbott et all I2OI0I ) and considering the case of (Tmax = 10 ^. The three lines correspond to different 
values of X = Rc/R, where Rc is the radius of the crystalline color superconducting core. 



which is defined in terms of the star principal moments of inertia. Standard neutron stars are expected to have 
maximum ellipticity of the order of Cmax ~ 10^^, meaning that this is the largest deformation sustainable by the 
crust before breaking. For the three-flavor crystalline phases considered in Section IlILDl deformed to the maximum 
extent allowed by the shear mo dulus, the maximum equatorial ellipticity sustainable could be as large as Cmax ~ 10"^ 



(IHaskell et all l2007t iLinl . 120071 ) . Different models of exotic stars composed of quark-clust ers have a m agnitude of the 



maximum ellipticity which is smaller than the CCSC phase by two orders of magnitude (jOwenl 120051) . 

The gravitational waves emitted by an optimally oriented compact star spinning at a frequency would produce 
a strain on Earth- based intereferometric detectors (such as LIGO, VIRGO, GEO600 and TAMA300) given by 



ho 



(200) 



where G is the gravitational constant and r is the distance. If the standard value Izz = 10^^ kg m^ is assumed, the 
above equation allows to relate the ellipticity to the spinning frequency and distance. The LIGO non-detection of 
gravity waves from nearby neutron stars already limits the pa, r amete r spa ce o f the m odel. A first analysis of CSOs 
which include the cr ystalline phase was made in Haskell et al\ (2007) andlLiij (2007) using the results of the S3/S4 
LIGO scientific run ( Abbott et al . 2007t ) and in Knippel and Sedrakian (|2009^ using the EoS discussed in Sec. IIV.DI 
and the results of the S5 LIGO scientific run ( Abbott et all . I2OO8I ) on gravitational wave emission from the Crab 



pulsar. In Fig. [T7]we present the exclusion plots for the Crab (left panel) and for the pulsar J2124-3358 (right panel), 
assuming maximum strain and that these stars are cons titut e d by incompressible and uniformly distributed matter. 
The dashed lines on the left panel are those obtained in iLfi] (l2007h for the Crab pulsar using the S3/S4 LIGO data 
assuming that the whole s tar is in the crystall ine color superconducting phase. We have determined the solid lines 
using the S5 LIGO data of Abbott et al. ( 2010t ). The parameter space results to be extremely reduced and even more 
severe constraints can be obtained from the pulsar J2124-3358, which leads to the exclusion plots reported on the 
right panel of Fig. [17] considering three different values oi x = Rc/R, where Rc is the CCSC core radius. 

Unfortunately, the maximum deformation of a star depends as well on the breaking strain, (Tmaxj which measures the 
largest shear stress deformation sustainable by a rigid body before breaking. For standard neutron stars it is assumed 
that 10~^ < (Tmax ^ 10~^. The results shown on the left panel of Fig. [T7] have been obtained with (Tmax = 10"'^ 
(upper red dashed line) and (Tmax = 10"^ (lower blue dashed line). The plots on the right panel are instead obtained 
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assuming (Tmax = 10~^. Considering lower values of (Tmax would move the curves up. It is important to remark that 
there seems to be little motivation for assuming that the breaking strain of the CCSC phase is of the same order 
of standard neutron star crust and only a study of the breaking mechanism of the crystalline pattern may allow to 
estimate this parameter. 

From Fig. [17] we infer that the S5 LIGO data give some hints about the properties of quark matter (under the 
assumption discussed above), because if it is present it should have a breaking strain (Tmax < 10^'^ and that if 
Cmax ~ 10"'^, it is unlikely that the pulsar J2124-3358 is a crystalline quark star; although, it could still contain a 
small crystalline quark core. With Advanced LIGO and Virgo detectors the spin-down limit on GW emission from 
known CSOs wi h reach e » 10 ~^, and objects with ellipticity of the order of 10~^ would be detectable up to the 
Galactic center (jPalombal 120121 ). These observation will put even more severe limitations on the parameters of the 
crystalline phases and on the EoS of highly compressed matter. 



B. Glitches 



The ri gidity of the crystal l ine phase may al s o be put in relation with the anomalies in the frequency of rotation 



of CSOs (jAlford et all 120011 : iMannarelli et a^.l . 120071) . A spinning neutron star observed as a pulsar gradually spins 



down as it loses rotational energy to electromagnetic radiation. But, every once in a while the angular velocity at the 
crust of the star is observed to increase suddenly in a dramatic event called a glitch. 



Pu lsar glitches are rare events fi rst observed in t he Vela pulsar ( Radhakrishnan and Manchester , 1969l : lReichlev and Downs 
1969) and in the Cra b pulsar ( Bovnton et all . Il969t Richards et all . 1969h which manifest in a variety of sizes 
(jEspinoza et ai . 2011), with an activity changing with time and peaking for pulsars with a characteristic age of about 
10 kyr. Although a numb er of models have been prop osed including two- fluid models (lAlpail. 1977t [Anderson and Itohl . 
1975 : Bavm et all. 1969) and "crust-quake" models ( Bavm and Pines . 1971 : Rudernianl . 1969l )(for a review of early 
models, see iRudermanI (jl972l )). these phenomena still remain to be well unde rstood , and their underlying mechanism 
is not yet completely clear. The standard explanation of pulsar glitches ( Alpai . 1977; Anderson and Itohl . 19751 ) 
requires the presence of two basic ingredients: a supcrfluid in some region of the star and a rigid structure that can 
pin the vortices without deforming when the vortices are under tension. As a spinning pulsar slowly loses angular 
momentum by radiation emission, the superfluid vortices tend to move toward the surface of the star — this happens 
because the angular momentum of a superfluid is proportional to the density of vortices — but, if the vortices are 
pinned to a rigid structure, they cannot move and after, as time passes, the lag between the superfluid component 
of the star and the rest of the star increases. This state persists until the stress exerted by the pinned vortices 
reaches a critical value, then there is a sudden avalanche in which vortices unpin, move outwards reducing the angular 
momentum of the supcrfluid, and then re-pin. As this superfluid suddenly loses angular momentum, the rest of the 
star, including in particular the surface whose angular velocity is observed, speeds up. 

The crystalline phases are rigid as well as supcrfluid, and it is reasonable to expect that the supcrfluid vortices 
within the rotating CCSC matter will have lower free energy if they are centered along the intersections of the nodal 
planes of the underlying crystal structure, i. e. along lines along which the condensate vanishes eve n in the absence of 



a rot ational vortex. A crude estimate of the pinning force on vortices within CCSC quark matter (jMannarelli et at. 



2007) indicates that it is comparable to that on neutron superfluid vortices within a conventional neutron star crust 
( Alpar et ai . 1984al lb). So, the basic requirements for superfluid vortices pinning to a rigid structure are all present. 
The central questions that remain to be addressed are the explicit construction of vortices in the crystalline phase 
and the calculation of their pinning force, as well as the calculation of the timescale over which sudden changes in the 
angular momentum of the core are communicated to the surface, presumably either via the common electron fluid 
or via magnetic stresses. Much theoretical work remains before the hypothesis that pulsar glitches originate within 
a CCSC core is developed fully enough to allow it to confront data on the magnitudes, relaxation timescales, and 
repeat rates of the observational data. 

The same mechanism outlined above might work for a star with a CFL inner core and an outer core in the CCSC 
phase. In this case superfluid vortices lying in the CFL phase are pinned to the periodic structure of the crystalline 
phase. 



C. Cooling and Urea processes 



Neutron stars cool down by neutrino emission and by photon emission, the latter dominating at late ages {t > 10^ 
yr). The cooling rate of a CSO may give information on its interior constitution, because different phases of hadronic 
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matter cool down in a rather different way, see e.g. iPetfiickl (|l992l ) for a brief review. 
The neutron star cooling is governed by the following differential equation: 

dT _ L, + _ Vne'l + V^el + L^ 
dt ^ 



V 



(201) 



where T is the inner temperature at time t and Ly and are respectively neutrino and photon luminosities, i.e. heat 
losses per unit time. After a very short epoch, when the temperature of the compac t star is of the order of ^ lO^^K 



and neutrinos are trappe d in the stellar core (jPrakash et ai . I2OOII : iRuester et all l2006bl : IShapiro and Teukolskvi 



1983 : Steiner et al. . 2002h . the temperature drops and the mean free path of neutrinos becomes larger than the star 
radius. For this reason the neutrino luminosity is obtained multiplying the emissivity by the corresponding volume, 
however one must distinguish between the emissivity e" of nuclear matter, which is present in the volume Ki, from the 
emissivity due to quarks, ej} if they are present in some volume V^; Cy and Cy denote specific heats of the two forms 
of ha dronic matter. We as sume a c ommon inner te mperature T , which is appropriate for sufficiently old compact 
stars ( Lattimer et al ]. ll994h . see 6.17. IHo et al\ (|2012l ) for a recent discussion. 

The luminosity by photon emission is instead a surface effect and can be estimated by the black-body expression 



(202) 



where R is the radius of the star, a is the Stefan-Boltzmann constant and the surface temperature is given by 

1/4 

T, -0.87 X 10^ ( I K, (203) 
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see 



/ rrn \ 0.55 

(ilk) 

T is the temperature at the basis of the stellar 



Gudmundsson et all (Il982h and IPage et~ai\ (|2004l ). where Ti, 
envelope and = GjqM/R'^ is the surface gravity. 

For what concerns the neutrino luminosity one has to consider the relevant weak processes. When kinematically 
allowed, direct Urea processes are the most efficient cooling mechanism for a CSO in the early stage of its lifetime 

However, the neutrino emission via the (nuclear ) direc t Urea processes, n ^ p + 
is only allowed for certain EoS (iLattimer et al\. 1991 ) having a sufficiently larg e 



(IShapiro and Teukolskv 
e + Vf._ and e~ + p 



1983) 



n -\- Vf, 



proton abundance to guarantee energy- momentum conservation ( Bahcall and Woh . Il965 : Chiu and Salpeteil . Il964 
Therefore, considering nuclear matter, only modified Urea processes are in general considered, where a bystander 
particle allows energy-momentum conservation. The resulting cooling is less rapid and the emissivity turns out to be 

2/3 

(204) 



(1.2 X lO'^crg cm-^s^i) — 



no 



10^ K 



see e.g. (jShapiro and Teukolskvll 19831) . much smaller than the emission rate e" — due to direct Urea processes. 
Here n is the number density and no = 0.16 fm^"^ is the nuclear equilibrium density. 

These considerati ons apply to stars containing standard nuclea r matter; faster cool i ng ca n be deter mined by the 



prese nce of a pion ([Bahcall and Wolfl . Il965: Maxwell et al\ . 119771 : iMuto and Tatsumi 119881 ) or kaon ([Brown et al. 



19881 ) condensate. Moreover, if the central region of the star consists of deconfined quark matter direct Urea processes 
involving quarks, i.e. the processes d ^ u + e~ -\- Ve and u + e~ ^ d -\- ma y take place and largely contribute to 



the cooling rate. It has been shown by Iwamoto ( Iwamotol 1980l 1981 . 1982) that quark direct Urea processes are 
kinematically allowed and the corresponding emission rate for massless quarks is of the order olsT^, where as is the 
strong coupling constant. This result is valid if quark matter is a normal Fermi liquid, but in the color superconducting 
phase the expression above is not correct because quarks form Cooper pairs and fermionic excitations are gapped. 
If the color superconductor is uniform, the corresponding neutrino emissivity and specific heat are suppressed by a 
factor e""^/"^, where A is the quasiparticle gap. This suppression is particularly strong in the CFL phase, because A 
is large and all quarks are paired. 

However, the ground state of quark matter in realistic conditions might not be the CFL phase, but the 2SC or the 
CCSC phase. In the following we shall evaluate the contribution of the latter phase to the neutrino emissivity and 
specific heat. 



1. Neutrino emissivity 

The transition rate for the /3-decay of a down quark da, of color a ^ r,g, b, into an up quark Ua 

daiPl) Ve{p2) + Ua{p3.) + e~ {pi) 



(205) 
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Wfi = V{2n)H\pi-p2-P3-p4)\M\^l[ 



1 



1=1 



2E,V 



(206) 



where V is the available volume and M is the invariant amplitude. Neglecting quark masses the squared invariant 
amplitude averaged over the initial spins and summed over spins in the final state is 



\M\^ = 64G|cos2 d,ipi-p2)iP3 ■ Pi) 



(207) 



where Gp is the Fermi constant and 6c the Cabibbo angle; we will neglect the strange-quark /3-decay whose contribution 
is smaller by a factor of tan^ 9c in comparison with (j207|) . Since for relatively aging stars there is no neutrino trapping, 
the neutrino momentum and energy are both of the order ksT. The magnitude of the other momenta is of the order 
of the corresponding Fermi momenta pi ^ pp ~ ^, ~ pp ~ fi and p4 ^ pp ^ ii^, which is smaller, but still sizable 



(see the discussion in Sec. IIII| . It follows that the momentum conservation can be implemented neglecting p2 and 
one can depict the 3-niomentum conservation for the decay ()205|) as a triangle ( Iwamoto . 1980l 1981 . 1982 ) having 
for sides pi , and P4 . It follows that we can approximate 



(Pi ■P2)(P3 •P4) ^ EiE2EzEi{l - cos6'i2)(l - cos 6*34) , 



(208) 



where Ej are the energies and O12 (resp. 6*34) is the angle between momenta of the down quark and the neutrino 
(resp. between the up quark and the electron). In the color superconducting phase one has to take into account that 
the neutrino emissivity 
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E2 M/fin(pi) [1 - n(p3)] [1 - n{pi)] {pi)Bl^ (p^) 



(209) 



depends on the Bog olyubov coefficients and B^^ which are functions of the quasiparticle dispersion laws 

(jAlford et aLl . l2005bl ). In Ea. (|209p the quark thermal equilibrium Fermi distributions 



n{pj) 



1 + exp 



T 



(210) 



appears because strong and electromagnetic processes establish thermal equilibrium much faster than weak interac- 
tions. The overall factor of 2 in Eq. (|209p keeps into account the electron capture proces s. 

The cooling of the CCSC matter with condensate (|150p was studied in Anglani et al. ( 2006f ). and in this case the 
largest contribution to the emissivity stems from the phase space region around the quark gapless modes, while the 
relevant momentum for the electron is its Fermi momentum, thus we have that 



d-'pi 



d^P3 



(TPA 



/^g dpi dVli Pf dpi dfliP^ dp3 dfl3 



(211) 



with dflj 



sin -dj ddj 



and Pi (resp. P3) is the quark down (resp. quark up) momentum where the corresponding 



quasi-particle energy vanishes. The gapless momenta Pi and P3 depend on the angle dj that quark momenta form 
with the pair momentum 2 q. In order to simplify the expression of the integral in Eq. (|209p we expand around the 
gapless modes, Ej{p) ~ fij + Vj [p — Pj) (for j = 1 and 3), with the quasiparticle velocity given by 



dE^ 
dp 



(212) 



In the three-flavor case the dispersion law of each quasiparticles has from one to three gapless modes, thus one has to 
expand the corresponding dispersion laws around each gapless momentum. 

Employing the above approximations the neutrino emissivity for each pair of gapless momenta Pi, P3, can be written 

as 



Gj, cos^ 0c fJ-l 
32^ 
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S'^^Hpi - P3-P4 -q){l- cos 012) (1 - cos 034) , (213) 



where I 
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Some of the angular integrations appearing in (j213p can be performed analytically, see lAnglani et al. 



(|2006f ) for more details, and the numerical computation can be reduced observing that even if each quasiparticle 
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dispersion law is characterized by various gapless momenta, not all of them satisfy the conservation of momentum 

Pi - P3 - P4 - g = 0. 

From the expression (|213p one can deduce that the largest contribution to the emissivity is due to blue quarks, 
that is to the process in Eq. (|205p with a = b. The reason is that according to the Ginzburg-Landau analysis of 



Casalbuoni et al\ (j2005al ) the gap parameter Ai vanishes and therefore the down blue quark is ungapped (see e.g. 
the discussion after Ea. (|12ip '). Then, the corresponding dispersion law is Ei{p) = p, with gapless momentum Pi ~ /i 
independent of mixing coefficient Bd^{Pi) = 1 and vi = \. The up-blue quark is instead paired with the stra, nge- 



red quark with gap parameter A2 = A and has two gapless momenta depending on 1^3, see Anglani et al. ( 20061 ) for 



an explicit expression. Of these two gapless momenta one satisfies the momentum conservation and thus the decay 
of a down-blue quark in an up-blue quark is not suppressed with respect to the analogous decay in ungapped quark 
matter, giving the leading contribution to the emissivity. Indeed, the case of quarks with colors red and green can be 
treated in a similar way, but for these colors neither the down nor the up quarks are un paired and the corres ponding 
processes are thus suppressed. These results are confirmed by the numerical analysis of Anglani et al. ( 2006f l. 



2. Specific heats 

As already discussed in Sec. III.DI for the two-flavor case, at low temperature the largest contribution to specific heat 
is determined by the fermionic quasi-particles. The specific heat of the three-flavor crystalline phase is given by the 
same formal expression given in Eq. (|67p. but with the quasiparticle dispersion laws of the three-flavor CCSC phase. 
The computation can be simplified following the same reasoning used above: the contributions of gapped modes are 
exponentially suppressed and each gapless mode contributes by a factor oc T. This result follows from the evaluation 
of the integral in Eq. (j67p employing the saddle point method and assuming that the quasi-particle dispersion laws 
are linear close the gapless momenta. Then, also the angular integral can be simplified because the dispersion laws 



arc gapless only in a restricted angular region. The numerical analysis confirms the above results, see lAnglani et al. 
([2OO6 ) for a discussion and for an expression of the specific heat. 

For unpaired nuclear matter and unpaired quark matter the contribution of each fermionic species can be approxi- 
mated by the fermionic ideal gas result 

klT 



^^PF^m^d^ + ippY , (214) 

where m and pp are the appropriate fermionic mass and Fermi momentum, respectively. For unpaired nuclear matter, 
the three sp ecies are neutro ns, prot ons an d electrons with Fermi momenta evaluated as in neutral matter in weak 



equilibrium ([Shapiro and Teukolskvt 119831 ) : 

1/3 / \ 2/3 



Pf 



(340 MeV) (^^^ , pI=p)^ (60 MeV) (^-^^ . (215) 

For neutral unpaired quark matter in weak equilibrium, the nine quark species have Fermi momenta independent of 
color and given by pp = /i H — , Pp = /i and pp = fi 



3. Cooling by neutrino emission 



In order to evaluate the effect on the compact star cooling of the CCSC phase we consider three different star toy 



models comprising neutral and /3 equilibrated matter ([Anglani et aLl . 120061 ). Model I is a star consisting of unpaired 



"nuclear" matter (neutrons, protons and electrons) with mass M = 1.4M0 (where Mq is the solar mass), radius 
i? = 12 km and uniform number density n = 1.5 no . Model II is a star containing a core of radius Rc = 5 km 
of unpaired quark matter with fi = 500 MeV, with a mantle of unpaired nuclear matter with uniform density n. 
Assuming a star mass M = 1.4 Mq from the solution of the Tolman-Oppenheimer-Volkov equations one gets a star 
radius i? = 10 km. Model III is a compact star containing a core of electric and color neutral three-flavor quark matter 
in the crystalline phase with gap parameter given in (|150p . with A2 = A3 ~ 6 MeV and fj, = 500 MeV, ml/^ = 140 
MeV. The outer part of the star is made of unpaired nuclear matter. Since the value of the gap parameter in the 
Model III is small, the radius of the CSO and of the quark core do not differ appreciably from those of a star with a 
core of unpaired quark matter, i.e. i?c = 5 km and i? = 10 km (also in these cases M = I.4M0). 

In Fig. [THI the cooling curves of the surface temperature as a function of time for the various models of star are 
shown. The solid black line refers to Model I; the dashed red line refers to Model II; the dotted blue line refers to 
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FIG. 18 (Color online). Surface temperature, in Kelvin, as a function of time, in years, for three toy models of pulsars. Solid 
black curve refers to model I; dashed red curve refers to model II; dotted curve blue curve refers to model III. Model I is a 
neutron star formed by nuclear matter with uniform density n = 0.24 fm~^ and radius R — 12 Km; model II corresponds to a 
star with R = 10 km, having a mantle of nuclear matter and a core of radius Rv ~ 5 Km of unpaired quark matter, interacting 
via gluon exchange; model III is like model II, but in the core there is quark matter in the LOFF sta te, see text for more d etails. 
All stars have M — IAMq. Parameters for the core are /i = 500 MeV and m?/^i= 140 MeV. From lAnglani et al\ (|2006D . 



Model III. Regarding the Model III, note that with increasing values of Mg the neutrino emissivity decreases (this 
is due to the fact that A decreases as one approaches the second order phase transition to the normal state) and 
in this case the quark matter in the star tends to become an unpaired Fermi liquid, for which the description of 
Iwamoto which includes Fermi liquid effects, should be used. For unpaired quark matter we use ~ 1, the value 
corresponding to fj. = 500 MeV and Aqcd = 250 MeV. The use of perturbative QCD at such small momentum scales 
is however questionable. Therefore the results for the Model II should be considered with some caution and the curve 
is plotted only to allow a comparison with the other models. In particular the similarity between the LOFF curve and 
the unpaired quark curve follows from the fact that in the CCSC phase the quasi-particle dispersion laws are linear 
and gapless, so that the scaling laws cy ^ T and e,y ^ arc analogous to those of the unpaired quark matter. 

From Fig. [15] we can see stars with a CCSC core cool down faster than ordinary neutron stars. This might 
have interesting phenomenological consequences because observational results on the cooling of pulsars are being 
accumulated at an increasing rate. Some data indicate that stars with an age in the range 10'^ — 10** years have a 
temperature significantly smaller than expected on the basis of the modified Urea processes. It is difficult however 
to infer, from these data, predictions on the star composition, as theses stars may have different masses. But, quite 
recently, the rapid cooling of the neutron star in Cassiopeia A ( Heinke a,nd Hoi 2010t Shternin et al . 201 ll) has been 
observed and explained as an effect of neutron superfluidity ( Page et all 2004 , 12011 ; IShternin et al. , l201l[) ; it would 
be interesting to check whether a similar behavior might be due to the presence of the CCSC phase. Note that the 
rapid cooling of neutron superfluidity is due to difermion pair-breaking effect, which enhances the emissivity of nuclear 
matter; the presence of blocking regions has a very similar effect for quark matter paired in a crystalline pattern. 

The analysis presented above is not conclusive because the thermodynamically favored phase should have a crys- 
talline pattern and not the simple modulation of (|150p . for which, however, the identification of the quasiparticle 
dispersion laws is still lacking. Nevertheless some qualitative assessments can be made from the obtained results. 
Slow cooling is typical of stars containing only nuclear matter or of stars with a uniform color superconducting phase 
(like CFL), thus the observed cooling rate of the neutron star in Cassiopeia A ( Heinke and Ho . 2010l : Shternin et al. 



20111) seems to disfavor both possibilities. This is compatible with the result that for intermediate densities the quark 



normal state and the CFL are less favored than the CCSC phase for a wide range of values of the baryonic chemical 
potential fj, and the strange quark mass Mg. Thus the presence of crystalline quark matter should be accompanied 
by a rapid cooling of the CSO. Note that the fast cooling of relatively young stars with a CCSC core is a consequence 
of the scaling laws for neutrino emissivity and specific heat. They depend on the existence of gapless points and 
follow from the existence of blocking regions in momentum space. Since this property is typical of any CCSC phase, 
independently of detailed form of the condensate, a rapid cooling should be appropriate not only for the simple ansatz 
assumed in Eq. (|150p . but also for condensates consisting of more plane waves. Note however, that the results reported 
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above do not take into account the transport properties inside the star, because it is assumed that at any time the 
external temperature is a given function of the internal one. It would be intere sting to include in the analysis the 



transport properties to have a detailed simulations of the CSO cooling, as e.g. in lHo et al\ (|2012l ) 



D. Mass-radius relation 



Since any phase transition leads to a softening of the equation of state ( Haensel 2003 : Lattimer and PrakashL 2001 ). 
it is in general expected that hyb rid stars — having a quark matter core and an envelope of baryonic matter — should 



have mass M < 1.7Mq, see e.g. (jAlford and Reddvl . 120031: iBubaUa et aLl . 120041: iMaieron et aU 120041 ) . For CSOs with 



larger masses, the deconfinement phase transition from baryonic to quark matter would reduce the central pressure 
to the point of instability towards black hole collapse. 

Some evidence for massive neutron stars with M ^ 2A/q has been inferred from various astro nomical observations : A 
compact star may exist in the LMXB (Low Mass X-ray Binary) 4U 1636-536 with 2.0±0.1Mq jBarret et a/.l.l2005l). A 



measurement on the pulsar PSR B1516+02B in the Globular Cluster M5 gave M = 2.08±0.19Mq dFreire et a?.l.l2008l) 



The m illisecond pulsar J1614-2230 has a mass (1.97±0.04)AfQ accurately measured by Shapiro delay (IDemorest et al.. . 
201(1 . Although these observations seem to disfavor the presence of quark matter (jLogoteta et all 120121 ). the details 



of a stability analysis depend on the theoretical model employed for the description of the hadronic phase, of decon- 
finement, and of th e CSC matter. We sh all show below that in the presence of the crystalline phase large masses can 
indeed be reached (llppolito et aLl . |2008| ). The drawback, as we shall see, is that if AI ~ 2Mq CSOs have a CCSC 



core, then ordinary M ~ 1.4M0 CSOs unlikely have a CCSC core. 

The description of quark matter relies on three different kind of models: the MIT bag model, the NJL models and 
the chromodielectric models. The results obtained within these three models may differ in a sizable way, and do as 
well depend on the detailed form of the model considered. 

Recent phenomenological studies of hybrid stars based o n the MIT bag model where carried out using a generic 
parameterization of the quark m atter EoS inlAlford et all ( 2005al ) and it was shown that hybrid stars may actually 
masquerade as neutron stars. In lAlford et all (|2005al ) non-perturbative QCD correc tions to the equati on of state of 
the Fermi gas were parameterized in a rather general way, with reasonable estimates (jFraga et aZ.1 . 120011 ). Taking into 
account these corrections, stable hybrid stars containing CFL quark matter may exist with maximum mass of about 
2Mq. 

Various studies of very massive hybrid stars within the three-flavor NJL model displayed a general instability 
towards collapse into a black hole (jBuballa et ali . |2004[ ). Stable stars featuring the 2SC phase we re obtained with 
typic al maximum masses M ~ 1.7M0 assuming reasonable values of the constituent quark masses ( Shovkovv et all . 



20031) or by replacing the hard NJL cut-offs by soft form factors with parameters fitted to a certain set of data 



(jBlaschke et all l2007t iGrigorian et al\. 120041). Heavier objects can be obtained if a repulsive vector interaction is 



introduced in the NJL Lagrangian lKiahZeiaI|,[2Q03), which makes the equation of state stiffer, but at the same 
time, reducing the amount of deconfined quark matter in the core of the star. 

Another source of uncertainty comes from th e nuclear equation of sta t e at h igh densities, which can be constructed 



starting from a number of different principles (jSedrakian . 



20071: IWebeij . Il999l) and a large number of EoS have been 



proposed for hybrid star configurations. In the analysis presented below, of the considered EoS, only the stiffcst were 
found to be admissible for phase equilibrium between nuclear and CCSC matter. 

The inclusion of hyperonic matter in the EoS of compact stars although reasonable is certainly troublesome, because 
the nucleon-hyperon and hyp eron-hyperon inter actions are not well known, even though some progress has been done 
mainly by lattice simulations (Aoki et al. . 2 pi2|). a nd introduce new parameters. Hyperons would certainly soften the 
equation of state (Glendenning and Kettned . 2000l ). making the comparison with the observed 2Mq CSOs problematic 
(jLattimer and Prakashll2007t ). 



1. Matching the equation of states 



The self-consistent computation of the strange quark mass given in llppolito et ali (j2007l ) consents to obtain the 
pressure as a function of the quark chemical potential fi. Thus, varying the quark chemical potential, the phase 
equilibrium between the confined and the color-superconducting phase can be constructed 

A normalization of the qua rk pressure in the NJL raodel is obtained by requirin g that the pressure vanishes at 
zero density and temperature ( Buballa and Oertel . 1999t Sandin and Blaschkel[2007l) . In the terminolo gy of the MIT 



bag model, this is equivalent to a subtraction of the bag constant from the thermodynamic potential (lAlford et 
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FIG. 19 (Color online). Left panel: Pressure versus density for the various models. The hues labele d as nuclear 1 (dashed blue) 
and n uclear 2 (dotted blue) refer to the EoS based on the Dirac-Bruckner-Hartree-Fock approach (|Webeii Il999l : IWeber et all 
|2007| ). Model A (heavy black) and of Al (medium-light red) have the nuclear (low-density) EoS corresponding to nuclear 2. 
The EoS of model B (light cyan) has the nuclear (low-density) EoS corresponding to nuclear 2. At the deconfinement phase 
transition there is a jump in the density at constant pressure. Right panel: Mass-radius diagram for non-rotating configurations 
including the the lower pulsar mass bound from J0737-3039 and the 2 solar masses upper bound (dashed horizontal line) from 
various astronomical observations. The solid h orizontal black line a pproximately corresponds to the onset of the crystalline 
color superconducting phase for model B. From ITppolito et al\ (|2008D . 



2005al ). Since the value of the bag constant is related to confinement, which is absent in the NJL model, it appears 



reasonable changing its value, and hence the normalization of the pressure. The simplest option is to consider the 
case of a constant shift in t he asy mptotic value of the pressure; alternatives include the use of form factors for the bag 



const ant (jGrigorian et all 120041 ) or the use of the Polyakov loop at nonzero density (jFukushimal . l2004t iRatti et 



20061) . 



Regarding the matching between the nuclear and quark matter E oS, it can be performed without changing the bag 
constant in the strong coupling limit, corresponding to Gjy/Gs = 1 ( Ruester et al. . 2005 : Sandin and Blaschkel 2007 ) 
(see Eq. ([T5|) and the discussion after Eq. (UHl))- The transition from the confined phase to quark matter happens 
when the baryo-chemical potential at which the pressures of the two phases are equal, meaning that the chemical 
potential curves P(/i) for these phases cross. This happens only in the strong coupling case, while at intermediate 
coupling and weak coupling, corresponding respectively to Gd/Gs ^ 0.75 and Gd/Gs S 0.7, the P{fJ.) curves for the 
equations of state of nuclear and quark matter do not cross; the models are thus incompatible, meaning that they 
cannot describe the desired transition between nuclear and quark matter. 

In the following we shall report on the analysis of llppolito et all (|2Q08( ). where the low-density equation of state of 
nuclear matter and the high-density equation of state of the CCSC matter are matched at an interface via the Maxwell 
construction. In the analysis the high-density regime is described by two EoS for crystalline color superconductivity 
which differ by the normalization of pressure at zero density (or, equivalcntly, the value of the bag constant). Of the 
various nuclear EoS considered only two are suitable to match with the CCSC equations of state. These are shown 
i n Fig. [191 as dashed blue and dot ted blue lines. The selected equations of state are the hardest two in the collection 
(jWebeiUlOQi IWeber et a^.l . l2007^ . 

Three different hybrid models are thus considered. For the models A and Al the nuclear (low density) EoS is the 
same and corresponds to the curve labeled as "nuclear 1" in Fig. [TOl For the model B, the nuclear (low density) EoS 
corresponds to the one labeled as "nuclear 2" in Fig. [THl The model Al has a high density EoS determined by the 
CCSC matter and normalized such that the pressure vanishes at zero-density. For the models A and B the high density 
behavior is determined by the same CCSC EoS used for model A, but the zero density pressure is shifted by an amount 
6p = 10 MeV/fm'^. This is equivalent to a variation of the bag constant whose value is related to confinement, which 
is unspecified within the NJL model, and is uncertain in general (Another possibility is to set dp = 0, but varying 
the value of the constituent masses of the light quarks in the fit of the parameters of the NJL model; f or small values 



of th e light quark masses the matching between quark and nuclear equations of state is facilitated (jBuballa et 
20041) ). 
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Note that because of the Maxwell construction of the deconfincnient phase transition, there is a jump in the density 
at constant pressure as illustrated in the left panel of Fig. [T9l 



2. Results 

Given the equation of state, the spherically symmetric solutions of Einstein' s equations fo r self-gravitating fluids 
are given by the well-known Tolman-Opennhcimer-Volkoff equations (iOppcnhei mer and Volkofi . 1939; Tolman, 193^ 

A generic feature of these solutions is the existence of a maximum mass for any equation of state; as the central 
density is increased beyond the value corresponding to the maximum mass, the star becomes unstable towards collapse 
to a black hole. One criterion for the stability of a sequence of configurations is the requirement that dM/dpc > 0, 
meaning that the mass of the star should be an increasing function of the central density. At the point of instability 
the fundamental (pulsation) modes become unstable. If stability is regained at higher central densities, the modes by 
which the stars become unstable towards the eventual collapse belong to higher-order harmonics. 

For configurations constructed from a purely nuclear equation of state the stable sequence extends up to a maximum 
mass of the order 2 Mq; the value of the maximum mass is large, since the chosen equations of state are rather stiff. 
The hybrid configurations branch off from the nuclear configurations when the central density reaches that of the 
deconfinement phase transition. The jump in the density at constant pressure causes a plateau of marginal stability 
beyond the point where the hybrid stars bifurcate. This is followed by an unstable branch [dM/dpc < 0). Most 
importantly, the stability is regained at larger central densities: a stable branch of hybrid stars emerges in the range 
of central densities 1.3 < pc/(10^^g cm^'^) < 2.5. Since the models A and B feature the same high-density quark 
matter, whereas the models A and Al the same nuclear equation of state, it is seen that the effect of having different 
nuclear EoS (the models A and B) at intermediate densities is substantial. At the same time, the small shift 5p by 
which the models A and Al differ does not influence the masses of stable hybrid stars, although it is necessary for 
matching of nuclear and quark EoS in the models A and B. It is evident that there will exist purely nuclear and hybrid 
configurations with different central densities but the same masses. T his is reminiscent of the situation encountered 
in non-superconducting hybrid stars ( Glendenning and Kettner . 2000() : the second branch of hybrid stars was called 



twin, since for each hybrid star there always exists a counterpart with the same mass composed entirely of nuclear 
matter. 

The right panel of Fig. [19] displays the astronomical bounds on the masses of CSOs (horizontal dashed lines) along 
with the results obtained with the previously discussed models. The upper bound corresponds to a mass of about 
M ^ 2Mq, while the lower bo und on the neutron star mass (1.249 ± 0.001)Mo is inferred from the millisecond binary 
J0737-3039 (iLvne et allbooS . Note that both the hybrid stars and their nuclear counterparts have masses and radii 



within these bounds. The hybrid configurations are more compact than their nuclear counterparts, i.e. they have 
smaller radii (an exception are those configurations which belong to the metastable branch). Contrary to the case 
of self-bound quark stars, whose radii could be much smaller than the radii of purely nuclear stars, the differences 
between the radii of hybrid and nuclear stars are tiny (< 1 km) and cannot be used to distinguish these two classes 
by means of current astronomical observations. 

The reported results differ from those described in lAlford et al. ( 2007 ) and Grigorian et al. { 2004 ) which have a 



narrow region of stable hybrid stars, branching off at the bifurcation point from the purely nuclear sequence. Indeed, 
the reported results cover a broad range of central densities, comparable to those covered by purely nuclear equations 
of state. 

For sequences constructed from models A and Al there is a range of masses an d radii that correspon ds to a family of 



stars that is separated from the stable nuclear sequence by an instability region (jlppolito et all\200&i \. The sequences 
corresponding to model B do not show such an instability region. 

The models A and Al are consistent with the bounds of M ^ 2Mq and for the model B these bounds correspond to 
the stable configuration with the largest mass. It should be noted that canonical 1.4Mq CSOs will be purely nuclear 
if they arc described by the three considered models. According with the presented analysis the CCSC matter appears 
only in sufficiently massive CSOs, with mass M > 1.7Mq for model B, corresponding to the solid horizontal line in 
the right panel of Fig|T9l 
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ACRONYMS AND SYMBOLS 

For clarity we report below a list of the most used 

BCC: Body centered cube 

BCS: Bardecn-Coopcr-Schrieffer 

BEC: Bosc-Einstein condensate 

CCSC: Crystalline color superconducting 

CFL: Color-flavor locked 

CSC: Color superconducting 

CSO: Compact stellar object 

EoS: Equation of state 

FCC: Face centered cube 

FF: Fulde-Ferrell 

GL: Ginzburg-Landau 

GW: Gravitational waves 

HDET: High Density Effective Theory 

LO: Larkin-Ochinnikov 

LOFF: Larkin-Ochinnikov- Fulde-Ferrell 

NGB: Nambu-Goldstone boson 

NJL: Nambu-Jona Lasinio 

QCD: Quantum chromodynamics 

QGP: Quark-gluon plasma 

2SC: Two flavor color superconducting 

gCFL: Gaplcss color-flavor locked 

g2SC: Gapless two flavor color superconducting 



acronyms and symbols. 

a, b: Adjoint color indices 
iTj,k: Flavor indices 
s,t: Spin indices 

Gd- Quark-quark coupling constant 
Gs'- Quark-antiquark coupling constant 
P^: Fermi momentum 
Ms'. Constituent strange quark mass 
Mq-. Solar mass 

a,l3,^: Fundamental color indices 
5^: Chemical potential difference 
A: Pairing gap 
ri: Free-energy 

Quark chemical potential 
He'. Electron chemical potential 
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